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A  method  for  calculating  solar  irradiance  with  cloud  and  other  information  is  des¬ 
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PREFACE 


This  is  the  final  report  of  research  conducted  for  the  Atmospheric 
Optics  Branch  of  the  Optical  Physics  Division,  Air  Force  Geophysics 
Laboratory.  The  research  was  conducted  under  Contract  F  19628  83  K  0040 
from  15  September  1983  through  14  December  1985  and  was  coordinated  with 
Dr.  Edmond  M.  Dewan,  Scientific  Program  Officer. 

This  report  elaborates  on  aspects  of  the  work  described  in  Quarterly 
Status  Reports  (021 082-1 -S  through  021 082-1 0-S )  and  in  the  first  Annual 
Report  (021082-1)  submitted  in  October  1984.  In  addition,  an  indexed 
compilation  of  193  abstracts  of  articles  and  reports  on  optical  turbu¬ 
lence  for  the  period  1970-1982  was  submitted  in  January,  1984.  Supple¬ 
menting  these  reports  of  progress  was  a  preliminary  version  of  a  stepwise 
procedure  for  computing  optical  turbulence  from  temperature  and  wind 
measurements  at  4  and  14  meters  that  was  submitted  in  September,  1984. 

^The  objective  of  the  research  was  to  develop  methods  for  estimating 
optical  turbulence  in  the  atmospheric  boundary  layer  with  conventional 
meteorological  variables.  There  are  many  such  variables,  considered 
here  to  be  those  observed  and  reported  hourly  from  weather  stations 
throughout  the  world.  An  additional  requirement  imposed  by  the  research, 
however,  is  that  the  variables  chosen,  compared  to  the  others,  also  have 
the  most  pronounced  effects  on  optical  turbulence  in  the  atmospheric 
boundary  layer. 

The  two  variables  that  meet  the  requirements  of  conventionality 
and  effectiveness  are  cloudiness  and  wind  speed.  Both  are  measured  or 
estimated  at  least  once  each  hour,  they  are  reported  in  standard  weather 


observations  and  together,  they  determine  the  intensity  of  optical 
turbulence  in  the  atmospheric  boundary  layer.  The  amounts,  types  and 
heights  of  clouds,  for  example,  through  their  effects  both  on  (1)  the 
amount  of  solar  radiation  available  to  heat  the  ground  in  daytime  and 
(2)  the  net  exchange  of  infrared  radiation  at  the  ground  surface  at 
night,  directly  affect  the  mean  vertical  temperature  gradient  and, 
therefore,  the  vertical  gradient  of  refractive  index  in  the  boundary 
layer.  Wind  speed  and  its  vertical  distribution  affect  not  only  the 
magnitude  of  the  mean  vertical  temperature  gradient  but,  in  concert  with 
it,  produce  the  turbulent  fluctuations  in  refractive  index  that  cause 
optical  turbulence.  Profiles  of  wind  and  temperature  in  the  boundary 
layer,  however,  are  not  routinely  available.  Although  the  estimates  of 
optical  turbulence  with  relationships  from  conventional  variables  are 
less  accurate  than  those  with  similarity  relationships  from  vertical 
profiles,  they  have  the  advantage  of  enabling  operational  estimates  to 
be  made  when  and  where  they  could  not  be  made  otherwise. 

The  research  is  primarily  analytical,  and  the  results  are  obtained 
from  an  analysis  of  optical  and  meteorological  measurements  made  by 
others  in  various  locations.  It  is  guided  by  knowledge  gained  in  simi¬ 
lar  work  here  several  years  ago.  ^Data  sets  consisting  of  boundary- 
layer  measurements  of  the  structure  function  parameter  for  refractive 
index,  (measured  optically  or  calculated  wit  a  various  methods  of 

''  »  *3 

estimating  the  temperature  structure  parameter  Cj),  and  meteorological 

variables,  including  solar  radiation  are  analyzed.  The  work  concentrates 

on  describing  optical  turbulence  only  for  sky  conditions  known  to  produce  / 

2  2 

the  largest  values  of  Cn(Cj)  in  the  boundary  layer,  both  in  daytime  and 


iv 


7^at  night,  in  that  only  conditions  with  a  cloudless  sky  or  with  scattered 
clouds  (half  or  less  of  the  sky  cloud-covered)  are  analyzed.  In  addition 
to  the  restriction  in  terms  of  cloudiness,  the  results  are  applicable 
mainly  to  optical  turbulence  conditions  over  land  surfaces  with  small 
amounts  of  soil  moisture., 
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1 .  INTRODUCTION 


The  propagation  of  optical  and  infrared  radiation  through  the 
atmosphere  is  often  seriously  limited  by  commonplace  atmospheric  charac¬ 
teristics.  In  addition  to  obvious  effects  of  absorption  and  scattering 
by  precipitation,  clouds  or  turbidity  constituents,  there  are  the  more 
subtle  clear-air  diffraction  and  refraction  effects  caused  both  by 
turbulent  fluctuations  of  atmospheric  density  and  by  average  vertical 
density  gradients.  To  an  observer  viewing  a  distant  object,  for  example, 
these  effects  may  appear  as  changes  in  the  object's  brightness,  position, 
size  and  color.  Brightness  fluctuations  are  termed  scintillation,  the 
motion  effect  is  called  shimmer  and  optical  turbulence  refers  to  the 
phenomenon  in  general.  In  the  atmospheric  boundary  layer,  optical  tur¬ 
bulence  is  usually  most  pronounced  in  cloudless  conditions  when,  except 
for  attenuation  by  atmospheric  turbidity  constituents,  optical  propagation 
and  the  transmission  of  visible  contrast  between  objects  and  their  back¬ 
ground  would  be  optimum.  In  cloudless  weather,  therefore,  this  phenomenon 
often  limits  the  performance  of  electro-optical  tracking  and  homing  sys¬ 
tems  such  as  high  energy  laser  systems,  long-range  reconnaissance  systems, 
laser  communications  systems  and  compensated  imaging  systems. 

An  Air  Force  need  for  a  method  to  characterize  optical  turbulence 
with  conventional  meteorological  variables  led  to  the  work  described 
herein.  The  objective  of  the  work  is  to  model  the  dependence  of  optical 
turbulence,  as  parameterized  by  the  refractive  index  or  temperature 

A  A 

structure  function  parameters  C^(Cy),  in  the  atmospheric  boundary  layer 
in  terms  of  standard  meteorological  and  terrain  information.  The  struc¬ 
ture  function  concept  is  treated  comprehensively  by  Tatarski  (1961),  and 


its  physical  meaning  in  relation  to  optical  turbulence  is  developed  and 
interpreted  by  Dewan  (1980).  The  temperature  structure  parameter  is 
defined  as 


Cj  =  [T(x)  -  T(x+r)]2  R"2/3  ,  (1) 


where  T  is  temperature,  x  and  r  are  position  vectors,  R  is  the  scalar 

magnitude  of  r  and  the  overbar  indicates  an  ensemble  average.  R  is 

a  constant  in  the  inertial  subrange  of  turbulence  eddy  sizes.  If 

2 

the  turbulence  is  homogeneous  and  isotropic,  CT  depends  not  on  r,  but 

2 

only  on  R.  Equivalently,  can  be  defined  as  the  amplitude  coefficient 
for  the  Variance  spectrum  of  temperature,  <j>j(k),  in  the  inertial  subrange 
as: 

*T(k)  «  0  25  C*  k_5/3  ,  (2) 


wnere  k  is  wavenumber  for  turbulent  fluctuations,  usually  calculated  as 

2iff/u,  where  f  is  frequency  and  u  is  horizontal  wind  speed. 

2 

The  refractive  index  structure  parameter  Cn  can  be  calculated  with 
(Kohsiek,  1985): 


C 


2 

n 


*•  2 


ATAQ  r 
TQ  TQ 


(3) 


where  Aj  and  Ag  are  given  by  Hill  fct.  al.  (1980),  T  is  temperature  (k)  and 
Q  the  absolute  humidity  (g  cm  ).  C-p  Cq  and  Cjq  are  the  structure  para¬ 
meters  for  temperature,  humidity  and  correlated  temperature-humidity 
fluctuations,  respectively,  and  can  be  calculated  with  measured  tempera¬ 
ture  and  humidity  fluctuations  (Kohsiek,  ’985) .  Because  the  present  work 
is  concerned  with  optical  wavelengths  and  a  continental  environment  only, 
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effects  of  humidity  fluctuations  are  assumed  to  be  small  compared  to 

effects  of  temperature  fluctuations,  especially  for  values  of  C^.  >_ 

10"\^m"^  (Kunkel  and  Walters,  1982);  (Wesely,  1976)  and  (Kunkel  et  al . , 

2 

1981).  For  this  condition  C  can  be  written  as: 

n 


rl  _  /  79  x  10'V  r2 
Ln  "  '  t2  ] 


(4) 


where  p  is  station  pressure  in  millibars  and  T  is  temperature  (K). 

2 

1.1  Cn  boundary  layer  model  variables 

The  atmospheric  density  fluctuations  responsible  for  optical 
turbulence  are  the  result  of  incomplete  turbulent  mixing  of  thermally 
stratified  layers,  a  condition  that  is  common  in  the  lower  atmosphere 
because  the  ground  usually  becomes  warmer  than  the  air  in  daytime 
(lapse  condition)  and  colder  at  night  (inversion  condition).  In  general, 
thermally  neutral  stability  in  the  first  several  meters  over  homogeneous 
surfaces  exists  only  in  very  cloudy  and  windy  conditions  and  for  very 
brief  periods  near  sunrise  and  sunset  when  the  vertical  temperature 
gradient  changes  sign  (Portman  et  al.,  1962).  As  a  result,  optical  tur¬ 
bulence  in  the  boundary  layer  varies  in  a  complex  manner  with  thermal 
stability,  but  it  is  at  a  minimum  in  thermally  neutral  (adiabatic) 
conditions  and  at  a  maximum,  dependent  on  wind  and  surface  roughness 

conditions,  for  maximum  vertical  temperature  gradients. 

2 

The  variables  affecting  Cn  in  the  boundary  layer  and  the  relation¬ 
ships  among  them  that  are  used  as  a  framework  to  model  the  variation 
2 

of  Cn  with  conventional  meteorological  and  surface  information  are  shown 
in  Fig.  1.  In  Fig.  1,  potential  temperature  fluctuations  (e')  are  the 
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Fig.  1.  Optical  turbulence  variables  and  relationships 


difference  between  the  actual  (e)  and  average  (s’)  potential  temperature 
and  can  be  expressed  as: 

e'  =  e  -  q  ■  f -| (de/dz ,w ' ,z)  *  f2(de/dz,du/dz,z,z0)  *  f3(de/dz,uz,z0).  (5) 

As  given  by  f  ^ ,  temperature  fluctuations  at  a  height  z  depend  on  the  mean 
temperature  gradient  (de/dz)  and  vertical  wind  speed  fluctuations  (w1). 

As  given  by  f2>  the  structure  of  the  wind  and  its  effects  on  w‘  at  a 
height  z  depend  on:  (1)  buoyancy  effects  (de/dz ) ,  (2)  vertical  wind 
shear  (du/dz),  and  (3)  the  transfer  of  momentum  by  the  drag  of  surface 
elements,  or  roughness  (zQ).  Finally,  as  given  by  f3,  assuming  (1)  that 
the  wind  speed  is  zero  at  the  roughness  scale  height  (zQ),  (2)  that  a 

- — — - — — . ..  ..  » n n «  n  rj >■>! rj, rj; nAnA  v,  irf'V. 
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wind  profile  (du/dz)  can  be  expressed  for  any  height  interval  (z-Zq), 

and  (3)  that  Zq  can  be  determined  from  the  wind  profile  when  de/dz  - 

l°C/100m,  the  structure  and  strength  of  turbulent  fluctuations  that 

2 

determine  optical  Cn  in  the  boundary  layer  can  be  considered  dependent  on 
the  mean  vertical  temperature  gradient,  the  average  wind  speed  at  a  given 
height  (uz),  and  the  surface  roughness,  Zg. 

Wind,  roughness,  and  height  can  be  regarded  as  secondary  factors 
affecting  both  the  temperature  profile  and  turbulence.  They  are  mecha¬ 
nical  effects  relatable  through  wind  shear,  and  can  be  combined  with 
buoyancy  effects  in  terms  of  the  vertical  temperature  gradient  as  the 
Richardson  number,  a  non-dimensional  ratio  that  is  a  measure  of  the 
relative  intensity  of  turbulence.  Portman  et  al.  (1962)  showed  that  the 
intensity  of  optical  turbulence  for  a  given  mean  temperature  gradient  at  2  m 
could  be  characterized  by  the  Richardson  number  (Ri).  They  showed  that 
for  inversion  conditions,  there  was  an  apparent  discontinuity  near  Ri  % 
+0.35  that  indicated  a  transition  from  turbulent  motion  to  a  more  wave¬ 
like  air  motion  for  larger  Ri. 

2 

The  relationships  diagrammed  in  Fig.  1  show  that  the  behavior  of  Cn 
is  traceable  to  effects  of  unique  combinations  of  cloudiness  and  wind 
speed  near  the  ground.  An  increase  in  cloudiness,  by  inhibiting  radia¬ 
tive  heating  or  cooling  of  a  surface,  decreases  the  vertical  temperature 
gradient  if  horizontal  advection  effects  are  minimal.  Similarly,  an 
increase  in  wind  speed  also  acts  to  decrease  the  vertical  temperature 

gradient  through  mixing  action.  Cloudiness  and  wind  speed,  therefore, 

2 

are  the  main  conventional  meteorological  variables  used  for  the  Cp 
model  development  discussed  below. 
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1.2  Relevant  boundary-layer  characteristics 

1.2.1  Unstable  stratification 

Most  boundary-layer  models  for  unstable  stratification 
partition  the  atmosphere's  first  1-3  km  according  to  the  physical 
processes  that  determine  mean  and  turbulent  properties  in  particular 
regions  (Webb,  1984;  Carson  and  Smith,  1974).  A  daytime  convective 
boundary  layer  (CBL)  is  usually  considered  to  be  that  region  within 
which  strong  vertical  mixing  occurs  as  a  result  of  net  upward  heat 
transfer  caused  by  solar  heating  of  the  surface.  It  is  usually  capped 
by  a  comparatively  non-turbulent  stable  layer.  Results  of  three- 
dimensional  numerical  calculations  (Deardorff,  1972,1974),  atmospheric 
measurements  (Kaimal  et  al.,  1976)  and  laboratory  experiments  (Deardorff 
et  al.,  1969),  show  that  the  two  most  important  variables  controlling 
mean  and  turbulent  properties  in  particular  regions  of  the  mixed  layer, 
or  CBL,  are  w*,  the  convective  velocity  scale  and  z-,  the  depth  of  the 
mixed  layer,  according  to: 

w*  =  [g/T  (we)oz.]1/3  ,  (6) 

in  which  g/T  is  a  buoyancy  parameter  (g  is  gravitational  acceleration, 
and  T  is  ambient  temperature),  (w?)Q  the  surface  kinematic  heat  flux 
and  z.j  the  height  of  the  mixed  layer.  The  large  convective  eddies 
scale  with  z^  and  the  turbulence  velocity  scale  is  proportional  to  w*. 
Both  w*  and  z.  undergo  diurnal  variations  that  start  at  zero  near  sun¬ 
rise  and  reach  midday  values  of  about  2  ms“^  and  1-3  km,  respectively. 

Several  definitions  of  z^  that  have  appeared  in  the  literature 
are  summarized  here  because  although  the  primary  maximum  of  C^(Cy) 
is  near  the  ground,  the  region  around  z^  is  a  region 
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with  a  secondary  maximum.  Scaling  theory,  furthermore,  indicates 
that  Cn(z)  and  Cy(z)  can  be  approximated  with  universal  functions 
of  z/z.  (Kukharets  and  Tsvang,  1980).  Kaimal  et  al .  (1976)  and  Caughey 
and  Palmer  (1979)  define  z^  as  the  height  of  the  inversion  base  capping 
the  mixed  layer  and  Wyngaard  and  LeMone  (1980)  define  it  as  the  minimum 
of  virtual  temperature  flux  (w?v).  As  pointed  out  by  Zhou  et  al.  (1985), 
in  most  cases  the  two  definitions  yield  comparable  values,  but  in  those 
cases  with  shallow  and  highly  baroclinic  boundary  layers  that  were 
analyzed  by  them,  the  inversion  base  had  undulating  motions  above  it 
and  turbulent  motion  below,  indicative  of  the  top  of  the  mixed  layer. 

The  height  of  minimum  temperature  flux,  however,  did  not  have  these 
characteristics  and  it  was  about  30%  higher. 

Kukharets  and  Tsvang  (1977)  and  Gamo  et  al.  (1976)  indepen¬ 
dently  found  that  a  reliable  and  physically  sound  experimental  determina¬ 
tion  of  z.  could  be  made  from  measurements  of  the  turbulent  energy  dissi¬ 
pation  rate,  e.  They  measured  e(z)  with  aircraft  and  found  that  in 
unstable  stratification  e(z)  showed  little  variation  with  height  until, 
at  a  height  usually  above  1  km,  it  decreased  2-3  orders  of  magnitude  in 
a  layer  about  200-400  m  thick.  The  subsequent  model  of  Kukharets  and 
Tsvang  (1980)  of  Cn(z)  was  based  on  this  finding. 

It  is  evident  from  the  definitions  of  z..  given  above  that 
detailed  measurements  of  various  meteorological  variables  are  required 
to  determine  z.  accurately.  Smith  (1977),  however,  through  an  analysis 
of  data  from  well-known  micrometeorological  field  experiments  in  England 
(Malvern  Experiment,  1976),  the  United  States  (The  Minnesota  1973  Experi¬ 
ment  and  the  Kansas  1968  Experiment),  and  Australia  (The  1967  Wangara 
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Experiment)  developed  a  practical  method  for  estimating  z^  with  conven¬ 
tional  variables.  The  variables  were  time  of  day,  month,  cloud  amount 
and  wind  speed  at  10  m.  His  method  consists  of  the  stepwise  procedure 
shown  in  the  nomogram  in  Fig.  2  ,  which  was  modified  slightly  for  the 
present  work.  As  discussed  below, reasonable  estimates  of  zi  were  ob¬ 
tainable  in  many  cases,  but  in  others,  large  disparities  between  estimated 
and  measured  values  were  observed. 

The  CBL  can  be  divided  into  three  regions  in  terms  of  zi : 

the  surface,  mixed,  and  interfacial  (entrainment)  layers.  The  surface 

layer  is  within  the  first  few  meters.  Its  mean  and  turbulent 

properties  are  determined  by  surface  fluxes  of  momentum,  sensible  heat 

and  latent  heat  and  hence  by  air-to-ground  (bulk)  differences  of  wind 

speed,  temperature,  and  vapor  pressure,  respectively.  In  the  surface 

layer,  the  assumption  that  the  fluxes  are  nearly  constant  with  height  has  been 

studied  experimentally  and  confirmed,  for  the  most  part,  for  sufficiently 

long  {%  1  hour)  averaging  times  (Haugen  et  al.,  1971).  Similarity  theory 

allows  relevant  micrometeorological  properties,  including  structure 

parameters,  in  the  surface  layer  to  be  represented  in  terms  of  height, 

appropriate  scaling  parameters,  and  the  Monin-Obukhov  (1954)  stability 

length  (Wyngaard  and  LeMone,  1980).  From  the  work  of  Wyngaard  et  al. 

o 

(1971),  for  the  free  convection  case,  close  to  the  ground  scales 

-2/3  -4/3 

according  to  z  ,  but  rapidly  approaches  z  '  . 

The  mixed  layer  extends  from  the  top  of  the  surface  layer 
to  near  the  base  of  the  capping  inversion,  0.1  z.  <  z  <  0.8  z...  Average 
wind  speed  and  potential  temperature  are  nearly  invariant  in  the  mixed 
layer.  Its  properties  are  determined  by  fluxes  in  the  surface  layer 


Fig.  2  A  nomogram  for  estimating  the  depth  of  the  boundary  layer 
in  the  absence  of  marked  advective  effects  or  basic  changes 
in  weather  conditions.  The  marked  example  shows  how  the 
diagram  is  to  be  used.  Adapted  from  Smith  (1977). 
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and  entrainment  processes  in  the  interfacial  layer  above.  In  its  lower 
2  -4/3 

half,  Cj  scales  according  to  z  ,  Obukhov  (1960).  In  a  transition 

layer  beginning  near  0.5  z^,  this  scaling  begins  to  break  down  at  heights 
dependent  mainly  on  the  value  of  z..  because  of  effects  of  turbulent  heat 
transfer  across  the  interfacial  layer  (Fairall  et  al.,  1982). 

The  interfacial  layer  at  the  top  of  the  mixed  layer  (0.8  z^ 

<  z  <  1.2  z.j)  separates  the  turbulent  boundary  layer  from  the  less  turbu¬ 
lent,  thermally  stable  air  of  the  free  atmosphere.  Turbulent  processes 
at  work  in  the  interfacial  layer  erode  it  by  entrainment  at  a  rate  that 
depends  on  the  intensity  of  turbulence  and  the  work  it  must  do  to  erode 
the  layer.  The  former  depends  on  surface  fluxes  and  w*,  and  the  latter 
depends  primarily  on  the  increase  in  temperature  through  the  layer  and 
the  lapse  rate  above  it.  Having  a  direct  effect  on  entrainment  processes 
are  thermal  plumes  or  updrafts  that  are  surrounded  by  larger  downdraft 
regions,  both  of  which  extend  from  near  the  ground  into  the  interfacial 
layer  (Webb,  1984).  Because  of  the  mixing  of  the  warm  dry  air  from  above 
the  inversion  with  cooler  and  moister  air  from  below,  secondary  maxima 
of  Cy  and  are  usually  observed  within  the  interfacial  layer  (Kukharets 
and  Tsvang,  1980;  and  Kaimal  et  al.,  1976). 

1.2.2  Stable  stratification 

Thermally  stable  conditions  in  the  boundary  layer  are  usually 
caused  by  nocturnal  infrared  radiation  loss  from  the  ground  and  atmos¬ 
phere,  if  the  sky  is  relatively  cloud-free,  but  they  can  also  be  caused 
by  advection  of  air  that  is  warmer  than  the  surface,  regardless  of 
cloudiness  or  time  of  day.  The  advection  of  air  above  freezing  over  a 
snow  or  ice  surface  is  an  example  of  the  latter.  In  contrast  to  the 
combined  effects  of  positive  buoyancy  and  vertical  wind  shear  that 
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enhance  vertical  motions  and  optical  turbulence  in  a  CBL,  there  are 
negative  buoyancy  effects  in  a  stable  boundary  layer  (SBL)  that  inhibit 
vertical  motion  and  reduce  the  length  scale  of  turbulent  fluctuations. 
When  the  length  scale  becomes  much  smaller  than  the  height  above  the 
surface,  effects  of  the  ground  on  turbulent  motion  decrease  to  the  point 
that  there  is  no  longer  an  explicit  dependence  on  height  (Nieuwstadt, 
1984).  Wyngaard  (1973)  expresses  the  behavior  of  surface-layer  profiles 
for  stable  conditions  in  terms  of  local  z-less  stratification. 

Results  of  field  measurements  made  here  (Ryznar  et  al.,  1971) 
and  observations  reported  by  Garrat  (1982)  are  among  those  that  show 
that  a  SBL  becomes  quite  well  established  within  about  2  hours  after  a 
transition  from  upward  to  downward  net  heat  flux.  Quite  often,  a  layer 
of  maximum  wind  speed,  ccnwonly  called  the  nocturnal  low-level  jet, 
develops  near  the  top  of  the  SBL  (Bonner,  1968).  As  discussed  by 
Blackadar  (1957)  and  Thorpe  and  Guymer  (1977),  the  jet  represents  an 
overshoot  in  the  increase  of  wind  speed  after  the  daytime  drag  by  tur¬ 
bulent  coupling  with  the  ground  surface  is  removed. 

As  Webb  (1984)  points  out,  there  is  no  unique  way  of  defining  the 
top  of  the  SBL,  denoted  here  as  zn<  In  numerical  models  of  the  SBL 
zn  is  usually  defined  either  in  terms  of  a  flux  criterion  as  the  height 
at  which  the  heat  or  momentum  flux  decreases  to  a  small  fraction  (\  5^) 
of  its  surface  value  (Brost  and  Wyngaard,  1978)  or  as  the  height  at 
which  a  critical  gradient  Richardson  number,  usually  between  0.2  and 
0.3,  is  reached.  In  terms  of  a  temperature  profile,  zp  can  be  defined 
by  the  height  to  which  a  ground-based  inversion  extends,  above  which 
the  potential  temperature  changes  little  with  height.  In  comparing 
these  methods  of  defining  zn>  Nieuwstadt  and  Dricdonk$(1979)  point  out 


that  a  turbulence-limit  height  for  z^,  if  defined  by  Ri  ^  0.2,  for 
example,  tends  to  level  off,  but  a  temperature-related  height  usually 
continues  to  increase  to  well  above  the  jet  maximum  because  of  radia- 
tional  cooling. 

In  later  work,  Nieuwstadt  (1984a,  b)  develops  local  scaling 
methods  to  describe  profiles  and  turbulence  structure  above  the  surface 
layer,  where  the  constant  flux  assumption  and  similarity  theory  lose 
their  validity.  Local  scaling  involves  the  local  kinematic  stress,  t» 
and  the  local  temperature  flux  wa  rather  than  the  surface  layer  variables 
u*  and  a*  contained  in  the  Monin-Obukhov  similarity  theory.  With  measure¬ 
ments  of  obtained  with  an  acoustic  sounder,  furthermore,  Nieuwstadt 
(1984)  found  that  zn  on  cloudless  nights  could  be  estimated  well  with 
a  diagnostic  relationship  involving  only  measurements  of  the  wind  speed 
at  10  m  as: 


7  -  oa  .,3/2 
Zn  '  28  u10 


(7) 


The  coefficient  is  probably  representative  only  for  the  re¬ 
gion  in  which  the  measurements  were  made,  but  as  discussed  below,  this 
relationship  seems  to  have  more  general  applicability.  It  is  a  practical 
means  of  obtaining  information  on  in  cloudless  nocturnal  conditions 
where  only  wind  speed  at  10  m  is  measured. 
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2  2 

1.3  Observations  and  models  for  C  (Cf)  in  the  boundary  layer 

Work  by  Wyngaard  et  al  (1971);  Kaimal  et  al .  (1976);  Caughey  et  al „ 

(1979);  Wyngaard  and  LeMone  (1980);  Burk  (1980);  Davidson  et  al .  (1981); 

Walters  and  Kunkel  (1981);  Fairall  et  al .  (1982);  Kunkel  (1982);  Kunkel 

and  Walters  (1982  and  1983)  and  Kohsiek  (1985)  shows  that  in  general,  the 

2  2 

height  variation  of  Cy  or  Cn  in  the  CBL  fits  within  the  framework  of 

Honin-Obukhov  similarity  theory.  A  z~^  near-ground  dependence  and 
-4/3 

z  mixed  layer  dependence  are  corroborated  by  experimental  results,  but 
both  exponents  have  been  observed  for  z  <  2m,  depending  largely 
on  vertical  wind  near.  Near  z^,  however,  uncertainties  regarding  entrain¬ 
ment  energy  (-.oyancy,  surface  wind  shear,  inversion  wind  shear  and  waves 
on  the  iversion  layer  itself)  and  interface  thickness  can  cause  large 
errors  when  model  results  are  compared  with  experimental  results  (Fairall 
at  al.,  1932). 

Davidson  et  al .  (1981)  and  Kunkei  (1982)  had  success  in  describing 

2 

Cn  at  one  height  in  relation  to  external  factors  such  as  surface-to-air 
temperature  difference  and  wind  speed.  Kunkel 's  experimental  work  con¬ 
siders  both  terrain  and  stability  effects  and  in  a  later  work,  Kunkel 

p 

and  Walters  (1983)  modelled  the  diurnal  dependence  of  Cn  by  solving  an 

extremely  detailed  energy  balance  equation  for  ground  surface  temperature 

for  a  clear  sky,  dry  soil  condition.  The  solution  yields  fluxes  of 

2 

sensible  and  latent  heat  with  which  C  is  calculated  near  the  ground. 

The  model  simulates  the  processes  quite  well  for  both  thermally  unstable 
and  stable  conditions  and  is  designed  for  surfaces  with  little  or  no 
vegetation.  Similarly,  Wesely  and  Alcaraz  (1973)  calculated  Cn  with 
estimates  of  sensible  and  latent  heat  flux  components  of  the  surface 
energy  balance. 
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Many  experimental  and  theoretical  results  of  investigations  of  the 
2 

height  variation  of  Cp  are  presented  and  discussed  in  detail  in  Gurvich 

et  al.  (1976).  Fig.  3  is  taken  from  their  book  and  shows  a  summary  of 

these  results,  with  our  translation  of  the  legend  given  below  it.  In  Fig.  3, 
2 

values  of  Cp  are  on  the  ordinate  and  height  in  meters  is  on  the  abscissa. 

2 

Roman  numerals  I -IV  on  the  Cp  axis  are  categories  based  on  the  findings 

p 

of  Voyt  et  al .  (1973),  who  measured  Cy  and  temperature  with  aircraft  at 
heights  from  50  m  to  5000  m  over  various  types  of  terrain  in  conditions 

from  very  unstable  to  near  adiabatic  stratification.  They  found  that 

2 

Cy(z)  in  the  50-ni  tc  600-m  layer  was  determinable  more  from  the  magnitude 
2 

of  Cy  at  50  m  than  from  the  temperature  lapse  rate.  For  example,  in 

2 

terms  of  their  3  categories  of  Cy(50),  they  found  tha-: 

I  tor  C2(50)  >  100xl0‘6(°C)2cm"2/3,  C*(z)  %  z'4/3; 

II  for  C2(50)  \  10-80xl0"6(°C)2cm'2/3,  C2(z)  %  z"2/3 
up  to  about  500-800  nr,  and 

III  for  Cy ( 50)  \  10  k(0C)2cm"2^3,  C2(z)  %  constant. 


Their  calculations  of  the  total  atmospheric  modulation  transfer 


function  for  the  50-5000  m  layer  showed  that  for  category  I,  the  contri¬ 
bution  of  the  5 On- 100m  layer  alone  was  77%  of  the  total,  for  II  it  was 


220  and  for  III  it  was  25. 


2 

Two  models  that  describe  C  (z)  in  a  CBL  with  z-  scaling  are  those 
of  Kaimal  et  al .  (1976)  and  Kukharets  and  Tsvang  (1980).  The  Kaimal 
model  was  developed  in  terms  of  surface  heat  flux  but  Walters  and 
Kunkel  (1981)  modified  it  for  more  practical  use  by  establishing 
experimentally  that  for  sufficiently  unstable  stratification, 


Fig.  3  Experimental  data  on  the  vertical  profile  of  the 

structural  characteristics  of  the  index  of  refraction 
in  the  atmosphere  (taken  from  Gurvich  et  al.,  1976). 

а,  b  and  c  -  outer  curves  corresponding  to  groups  I,  II  and  III  of  the  data 

of  Voyt  et  al . ,  (1973) ; 

d  -  neutral  data  (Kail istratova ,  1966)  obtained  from  simultaneous  surface 
(z  =  4m)  and  aircraft  (z  =  50m)  measurements; 

e  -  measurements  with  resistance  thermometers  suspended  from  a  tethered 
balloon  (Coulman,  1969); 

f  -  data  obtained  on  a  meteorological  tower  (Byzova  and  Vyal’tseva,  1970) 
from  micropulsation  measurements  of  temperature; 

g  -  data  from  Bufton  et  al.  (1972)  opined  from  measurements  of  micropuisa 
tions  of  temperature  from  a  balloon  in  free  flight; 

h  -  estimate  of  from  astronomical  observations  (Tatarski,  1967); 

i  -  aircraft  optical  measurements  (Morris,  1973); 

1  -  C2(z)  w-4/3; 

2,  •>,  4  -  cjjlz)  ^  z'2/3i 

5  -  calculated  from  the  dissipation  of  turbulent  kinetic  energy  (Hufnagel , 

I  0* ?  ,  and 

б,  7  -  c‘  exp  [-  (z-z2)/z3]  . 
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a  knowledge  of  at  one  height  near  the  ground  and  information  on  2^ 


are  sufficient  for  describing  C^(z)  in  the  following  regions: 


cn<*) 


.{ 


(2/Zl) 


-4/3 


for  2^ ,  z  <_0.5  z. 


(8) 


,-4/3 


for  0.5  z.  <  z  <  0.7  ,  and 


(0.5  z./ Zj) 

2.9  (0.5  zi/z1)"‘t/J(z/z  .y  for  0.7  zi  <z<z. 


i-4/3/_ x3 


It  can  be  noted  that  the  validity  of  the  Kaimal  model  extends  only 

p 

as  high  as  and  in  terms  of  C  ,  the  model  shows  an  increase  with  height 
near  Zy 

The  model  of  Kukharets  and  Tsvang  (1980), hereafter  referred  to  as 

K-T,  is  based  on  an  analysis  of  2  years  (1976  and  1977)  of  measurements 
2 

of  CT(z)  and  meteorological  variables  in  the  CBL  over  steppe  and  forested 

2 

surfaces  in  the  Soviet  Union.  They  normalized  their  measurements  of  Cy 
2 

to  Cy(0.1),  which  corresponds  to  a  dimensionless  height  c  =  z/z^  =  0.1, 
where  z^  was  determined  from  the  abrupt  decrease  in  the  turbulence  energy 
dissipation  rate.  An  empirical  equation  was  developed  to  fit  a  composite 
of  observed  variations  of  a  dimensionless  structure  parameter  <i j  = 
Cy(;)/Cy(0.1 )  through  a  secondary  maximum  near  z/z^  %  1.1  without  a 
height-dependent  breakdown  as  in  the  Kaimal  model.  Their  equation  is 


.4,,  -Mt-l-l)2 

*  2U)  =  k-j  C  +  t^exp  6 


(9) 


with  the  values  of  the  coefficients  k^  =  4.6x10  ,  k2  =  0.6  and  k^  =  12. 

-4/3 

It  can  be  noted  that  eq.  9  describes  a  z  variation  at  altitudes  well 

p 

below  2^ ,  with  a  Gaussian-shaped  variation  of  C*  with  height  through  z^ . 
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The  coefficients  k^  and  k^»  respectively,  determine  the  amplitude  and 
standard  deviation  of  the  Gaussian-type  variation  through  z.. .  The 
larger  the  magnitude  of  k^,  for  example,  the  smaller  the  standard  devia¬ 
tion  of  the  Gaussian-shaped  variation. 

2  2 

Profiles  of  C  (CT)  calculated  with  the  modified  Kaimal  model  and 
n  T 

the  K-T  model  were  compared  with  balloon-borne  and  scintillometer 
measurements  of  profiles  for  conditions  of  unstable  stratification  by 
Murphy,  Dewan  and  Sheldon  (1985).  They  adapted  the  K-T  model  for  appli- 

p 

cation  to  measurements  of  Cn  at  a  reference  height  of  14  m.  With  profile 
values  of  z^ ,  they  found  reasonable  agreement  for  an  average  of  several 
profiles,  but  for  single  profiles,  they  found  that  measured  Cn(z)  in  the 
vicinity  of  z^  was  narrower  than  the  K-T  model  prediction.  They  suggested 
that  a  value  of  k^  =  580  instead  of  k^  =  12  is  more  applicable  to  single 
profiles. 

For  stable  stratification,  fewer  research  results  are  available  and 

less  modeling  work  has  been  conducted.  Caughey  et  al .  (1979)  reported 
2 

measurements  of  ay  for  heights  from  4  m  to  305  m  for  a  site  in  Minnesota. 

The  behavior  of  the  corresponding  structure  parameter  normalized  by 

the  surface  value  of  a  scaling  temperature  e*  as  Cy  zn'  /e*  decreased 

linearly  from  about  140  for  small  z/zn  to  near  unity  near  the  top  of 

SBL  (Caughey,  1982).  Their  results  also  indicate  that  for  stable 
2 

conditions,  Cy  can  be  described  in  terms  of  Monin-Obukhov  similarity 

in  the  lowest  few  meters  of  the  boundary  layer  if  it  is  assumed  that 

the  atmosphere  is  fully  turbulent.  Wyngaard  et  al.  (1971)  found  that  Cy 

decreases  with  height  more  slowly  than  z"2^3  for  stable  conditions  and 
-2/3 

approached  z  for  near  neutral  conditions. 
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2.  Data  acquisition,  processing  and  tabulation 

Early  in  the  work,  a  search  was  conducted  for  experimental  data 
that  included  simultaneous  measurements  of  structure  parameters  and 
meteorological  variables.  A  copy  of  the  PRESSURS  data  base  survey 
(AFWAL-TR-83 -1 077 )  conducted  for  the  EO  Sensor/ Atmospheric  Science  Group 
at  the  Air  Force  Avionics  Laboratory  was  obtained  that  led  to  the  iden¬ 
tification  of  data  sets  considered  to  be  complete  enough  for  analysis 
and  obtainable  in  a  format  compatible  with  computing  facilities  here. 

The  sources,  responsible  individuals  and  types  of  data  are  shown  in 
Table  1. 


Table  1  Structure  parameter  and  meteorological 
data:  sources,  types  and  responsible 
individuals 


Source 

1)  Air  Force  Geophysics  Laboratory 
E.M.  Dewan,  J.  Brown 

See  Brown  et  al .  (1982)  and 
Brown  and  Good  (1984). 


2)  Rome  Air  Development  Center 
Advanced  Optical  Test  Facility 
J.  Foster 


3)  NOAA/ERL  Wave  Propagation 
Laboratory  (Boulder  Atmosphe¬ 
ric  Observatory) 

J.  Gaynor 


Type  and  Form 

O 

a)  Balloon-borne  Cn  (thermosonde) 
and  met  variables  through  tro¬ 
posphere  for  several  locations 
(mag  tape). 

9 

b)  Optical  and  spaced-probe  C^,  mi- 
crometeorological  profiles  and 
radiation  variables  for  field 
experiments  at  White  Sands,  NM 
(graphs  of  diurnal  variations). 

a)  Optical  Cn  at  2-m  high,  118-m 
long  optical  path. 

b)  Met  variables  and  solar  radiation 
at  2  m  (3-min  averages  on  mag 
tape) . 

2 

a)  Optical  Cn  at  4  m  on  3  legs  of 
triangle,  450-m  optical  path. 

b)  Temperature  spectral  densities 
and  met  variables  for  10,  22,  50, 
100,  150,  200,  250  and  300  m  on 
300-m  BA0  tower.  Also  solar 
radiation,  (20-min  averages  on 
mag  tape) . 


2.1  Air  Force  Geophysics  Laboratory  (AFGL)  data 

Data  for  the  research  were  provided  by  AFGL  in  the  following  forms: 

(1)  computer  tapes  with  data  for  17  daytime  and  19  nighttime  soundings 
o 

of  C^,  temperature,  relative  humidity  and  wind  for  Westford,  MA,  Hanscom 
AFB,  MA;  and  Boulder,  CO  and  (2)  similar  data  plus  graphs  of  optical 
and  meteorological  variables  at  White  Sands,  NM.  Relevant  information 
was  obtained  from  (1),  but  because  of  the  greater  detail  and  completeness 
of  (2)  from  heights  near  the  ground  through  the  complete  boundary  layer, 
greater  emphasis  was  placed  on  an  analysis  and  interpretation  of  the 
latter  set  of  data. 

2 

The  meteorological  and  Cn  data  for  special  experiments  at  the  White 
Sands  Missile  Range  (WSMR)  in  New  Mexico  were  provided  for  the  research 
in  the  following  forms:  (1)  82  vertical  soundings  of  Cn,  temperature, 
relative  humidity  and  wind  on  magnetic  tape  (2)  graphs  of  diurnal 
variations  of  meteorological  and  optical  variables  in  the  atmospheric 
surface  layer  and  (3)  near-ground  measurements  of  optical  variables 
on  magnetic  tape.  The  vertical  soundings  had  20-m  height  resolution 
and  were  measured  with  balloon-borne  thermosondes  (Brown  et  al.,  1982) 
by  personnel  from  the  AFGL  Atmospheric  Optics  Branch,  Optical  Physics 
Division,  and  the  ground-based  measurements  were  made  by  personnel  from 
the  US  Army  Atmospheric  Sciences  Laboratory,  WSMR.  The  field  experiments 
for  which  data  were  obtained  were  entitled:  CRITICAL  LASER  ENHANCING 
ATMOSPHERIC  RESEARCH  (CLEAR),  and  took  place  from  27  Aug.  to  27  Sept. 

1984  (CLEAR  I),  from  25  Feb.  to  10  Mar.  1985  (CLEAR  II),  and  from  29  July 
to  3  Aug.  1985  (CLEAR  III). 

A  significant  effort  was  devoted  to  devising  a  method  for  programming 
the  University  of  Michigan  computer  to  accept  the  format  of  the  CLEAR  I 

- - - 
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and  II  sounding  data  for  processing  and  analysis.  Once  a  method  was 
developed,  the  Cn  data  for  each  sounding  were  averaged  logarithmically 
within  the  height-dependent  increments  shown  in  Table  2  .  This  method 
of  averaging  reduced  the  number  of  sample  points  in  3000  m  from  about 
150  to  23  and  distributed  them  quite  equally  on  a  logarithmic  scale. 

Table  3  shows  the  averaging  method  applied  to  the  temperature  and 
humidity  data,  which  were  averaged  linearly  in  the  height  intervals 
shown.  Most  boundary  layer  wind  measurements  were  missing. 

For  visually  evaluating  each  sounding  averaged  in  these  ways,  com- 

_2 

puter  graphs  were  produced  with  the  21  to  24  values  of  Cjj  plotted  against 
height  up  to  3000  m,  with  both  variables  in  logarithmic  coordinates,  and 
with  35  to  40  values  of  temperature  and  relative  humidity  up  to  3000  m, 
with  both  variables  in  linear  coordinates.  For  nocturnal  conditions, 
profiles  with  all  30  values  up  to  a  height  of  600  m  values  were  also 
graphed  to  provide  improved  resolution  in  determining  the  top  of  the 
stable  boundary  layer. 

Because  the  vertical  soundings  for  CLEAR  I  and  II  were  obtained  at 
various  times  of  day,  they  were  representative  of  various  conditions  of 
thermal  stability.  Stability  conditions  estimated  from  the  graphs  of 
diurnal  temperature  profiles  measured  on  HIDL  tower  2  at  heights  of  0.5, 
1,  2  and  4  m  were  used  to  group  only  those  soundings  obtained  in  steady 
cloudless  or  scattered  sky  conditions  into  general  categories  of  stable, 
unstable  and  near  neutral.  Table  4  contains  (1)  a  listing  of  the 
soundings,  (2)  the  lowest  height  of  the  first  measurement  (z^),  (3) 
cloud  amount  and  height  for  the  sounding  time  as  estimated  from  hourly 
weather  observations  at  C  Station,  WSMR,  and  at  Holloman  AFB,  (4)  wind 


o 

Table  2  Averaging  method  for  log  Cn  for  the  AFGL  data, 
Zj_  refers  to  the  lowest  measurement  height  on 
each  sounding 


Height  Interval  (m) 


Number  of  Values  Averaged 


zl-20 

1 

21-40 

1 

41  -60 

1 

61-80 

1 

81-100 

1 

101-120 

1 

121-140 

1 

141-160 

1 

161-200 

2 

201 -240 

2 

241 -280 

2 

281-320 

2 

321-380 

3 

381 -440 

3 

441-520 

4 

521 -640 

6 

641 -800 

8 

801-1000 

10 

1001-1240 

12 

1241-1560 

16 

1561-1960 

20 

1961-2440 

24 

2441-3000 

28 

wWy'vWJ 
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Table  3  Averaging  method  for  temperature  and  relative 
humidity  for  the  AFGL  data.  z\_  refers  to  the 
lowest  measurement  height  on  each  sounding 


Height  Interval  (m)  Number  of  Valties  Averaged 


zl-20 

21-40 

41-60 


1 

1 

1 


to  200  m  1 

201-240  2 

241-280  2 

281-320  2 

,  • 

•  * 

to  600  m  2 

601 -680  4 

681-760  4 

761-840  4 

•  • 

to  1400  m  4 

1401-1560  8 

1561-1720  8 

1721-1880  8 

•  • 

*  • 

to  3000  m  8 


2 

Table  4.  Values  of  meteorological  variables,  Cn  (scintillo¬ 
meter)  and  lowest  sounding  heights  for  thermosoundings 
of  CLEAR  I  and  CLEAR  II 


CLEAR  I 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

DAY/TIME 
1984  (MSI) 

STAB. 

2l 

M 

AMT 

CLOUD 

HI 

u(4m) 

ms"^ 

xlC 

C*(4m) 

,-14-2/3 
)  m 

SOLAR 

_2 
Wm  c 

zi 

(m) 

zn 

(m) 

REMARKS 

AUG  28/1354 

u 

12 

3 

60,120 

MISDA 

1400 

no  cl 

29/1502 

22 

9 

60,120 

2.6 

1.2 

145 

29/1753 

513 

9 

65,100 

2.6 

1.1 

50 

29/2104 

s 

20 

3 

80,250 

0.9 

7.5 

180 

31/1404 

u 

10 

3 

80,250 

2.9 

27 

745 

2600 

"oCn 

SEP  1/1525 

12 

7 

60,120 

1.4 

9 

500 

2/1909 

s 

22 

3 

50,100 

MISDA 

- 

COMPLEX  T 
PROFILE 

3/2113 

s 

16 

1 

100 

1.8 

4 

80 

4/1348 

u 

17 

1 

60 

1.6 

38 

830 

1800 

4/1716 

n 

13 

1 

200 

1.8 

2.5 

220 

5/1319 

u 

20 

0 

1.9 

51 

890 

1800 

5/1716 

n 

21 

0 

2.9 

2.1 

190 

6/1311 

u 

20 

0 

2.9 

44 

890 

1800 

6/2056 

s 

18 

1 

250 

1.0 

7 

no 

7/1417 

u 

13 

1 

250 

3.5 

34 

780 

2100 

7/1708 

n 

22 

0 

2.3 

4.4 

220 

7/1724 

n 

28 

0 

1.3 

2.4 

140 

7/1945 

s 

8 

0 

1.2 

18 

60 

8/1343 

u 

25 

2 

60,250 

1 .6 

46 

840 

2600 

8/2010 

s 

4 

1 

120,250 

1.0 

6.5 

140 

9/1303 

u 

15 

1 

120,250 

2.1 

20 

880 

1900 

9/1746 

n 

10 

1 

80,250 

0.6 

M 

10/1457 

u 

15 

3 

60,250 

3.1 

17 

570 

2700 

10/1733 

n 

1908 

4 

60,250 

1.5 

2.8 

95 

10/2005 

s 

12 

4 

60,250 

1.0 

13 

90 

no  cj; 

11/1303 

u 

37 

4 

65,200 

1.9 

47 

870 

2000 

11/1627 

u 

10 

3 

65,120 

2.6 

5.5 

225 

2300 

11/1854 

s 

8 

2 

65,100 

0.9 

4 

60 

12/1157 

u 

8 

1 

50,120 

1.5 

34 

885 

1000 

12/2016 

s 

18 

3 

60,120 

0.6 

7 

1 404 
1300 

DOUBLE  LAPSE 

(Continued  on  next  page.) 
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(Table  4,  continued.) 


DAY/TIME 
1984  (MST) 

O) 

STAB. 

(2) 

ZL 

M 

(3) 

CLOUD 

AMT  HT 

(4)  (5) 

u(4m)  Cn(4n,) 

ms’’  xlO-’V2'3 4 * 6 7 8 

(6) 

SOLAR 

_2 
Wm  c 

(7) 

zi 

i ml 

(8) 

zn 

M 

(9) 

REMARKS 

SEP  13/1301 

u 

n 

2 

70,100 

1.7 

46 

865 

2000 

13/1931 

s 

21 

3 

60,120 

1.3 

7 

190 

14/1258 

u 

26 

2 

60,120 

2.4 

44 

854 

2000 

14/1609 

u 

16 

2 

60,250 

3.1 

7.8 

411 

2700 

15/2019 

s 

17 

5 

60,120 

MISDA 

560 

16/1243 

u 

17 

1 

70 

2.5 

41 

750 

J  200 
’1500 

DOUBLE  zi 

16/1851 

s 

17 

2 

70,100 

1.7 

4.8 

60 

17/1854 

s 

18 

3 

50,250 

0.9 

3.7 

100 

18/1253 

u 

19 

3 

60,120 

1.8 

46 

750 

2600 

18/1851 

s 

7 

2 

60,120 

0.8 

7 

100 

20/0321 

s 

7 

0 

1.5 

5.8 

50 

20/1034 

u 

22 

3 

60,110 

2.2 

48 

810 

2800 

no  C* 

21/0302 

s 

22 

2 

100,200 

1.5 

9.1 

i  60 
<180 

DOUBLE  LAPSE 

21/0751 

n 

14 

3 

120,200 

1.0 

12 

380 

24/1239 

u 

22 

1 

60,250 

2.6 

19 

2400 

24/1908 

s 

8 

1 

60,250 

0.9 

1.7 

120 

25/1232 

77 

7 

60,120 

6.0 

22 

830 

26/0658 

930 

10 

12,20 

2.6 

0.02 

35 

26/1236 

39 

10 

8,22 

2.0 

0.8 

180 

27/1155 

3 

5 

40 

1.6 

15 

920 

(1)  Stability  category:  unstable  (u ) ,  neutral  (n),  stable  (s)  for  soundings  with 

cloud  amount  <  6/10. 

(2)  z^:  lowest  sounding  height. 

(3)  Cloud  amount  (tenths)  and  height  (hundreds  of  feet)  for  Holloman  AFB  (CLEAR  I). 

(4)  u(4m):  4-meter  wind  speed  for  Tower  I,  HIOL. 

(51  CM 4m):  4-meter  scintillometer  value,  Tower  I,  HIDL. 

(6)  SOLAR:  direct  plus  diffuse  irradiance  on  a  horizontal  surface,  Tower  I,  HIDL. 

(7)  z^:  top  of  mixed  layer. 

(8)  zn:  top  of  nocturnal  boundary  layer. 


Table  4  (Continued.) 

2 

CLEAR  II  (no  ground-based  meteorological  or  Cn  data  available.) 


(1) 

(2) 

(3) 

(7) 

(8) 

(9) 

DAY/TIME 

ZL 

CLOUD 

zi 

zn 

1935  (MST) 

STAB. 

ini 

AMT 

HT 

iml 

M. 

REMARKS 

FEB  26/1429 

9 

10 

280 

26/1759 

201 

8 

160,260 

I  450 
'1000 

27/0938 

u 

12 

3 

220 

DOUBLE  INVERSION 

27/1538 

12 

7 

80,220 

27/1805 

13 

10 

110,220 

27/2106 

s 

4 

4 

120 

400 

28/1403 

3 

6 

120,220 

28/1815 

s 

12 

2 

60,250 

70 

28/2104 

s 

10 

1 

60 

80 

MAR  1/0950 

u 

13 

1 

60 

1700 

1/1527 

n 

0 

2 

65,250 

j  40 
'160 

1/1807 

s 

14 

4 

65,240 

DOUBLE  LAPSE 

1/2059 

s 

14 

5 

250 

J  50 
'230 

DOUBLE  LAPSE 

2/0856 

7 

10 

150,250 

2/1736 

3 

8 

80,250 

4/0904 

n 

35 

CLR 

4/1756 

s 

15 

II 

I 

4/2050 

s 

4 

II 

70 

5/0901 

n 

1 

II 

5/1441 

u 

44 

1 

250 

1100 

HIGH  zL 

5/1746 

s 

159 

3 

150,250 

- 

HIGH  2l 

6/0850 

9 

10 

120,250 

6/1731 

204 

10 

150,200 

6/1959 

31 

9 

150,250 

7/0846 

31 

10 

150,200 

7/1452 

16 

10 

120,200 

7/1735 

3783 

10 

150,200 

7/2002 

6 

9 

150,220 

8/1311 

8 

8 

120,200 

8/1721 

13 

10 

150,250 

9/0854 

17 

10 

100,200 

9/1928 

s 

8 

5 

120 

- 

COMPLEX  T  PROFILE 
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p 

speed,  and  (5)  scintillometer  measurements  of  Cn-  The  values  of  wind 
2 

speed  and  CR  shown  in  Trble  4  are  those  measured  at  4  m  on  HIDL  Tower  1. 

Because  the  height  of  first  measurement  varied  with  each  sounding,  values 
2 

of  Cn  obtained  with  the  thermosonde  were  not  available  for  a  standard 
height.  Such  a  measurement  would  have  been  useful  in  permitting  a 

grouping  of  the  soundings  not  only  according  to  stability  criteria  but 

2 

also  according  to  values  of  Cn  at  one  height. 

Also  shown  in  Table  4  are  estimates  of  z.,  the  top  of  the  mixed 
layer  for  unstable  stratification,  and  zn,  the  top  of  the  nocturnal 
boundary  layer  for  stable  conditions.  Values  of  z.  and  zn  were  deter¬ 
mined  by  inspecting  qraphs  of  profiles  for  height  intervals  with  changes 
,  2 

in  temperature,  relative  humidity  and  Cn*  For  many  profiles,  z^  was 

detectable  as  a  temperature  inversion  above  an  adiabatic  layer,  but  a 

more  reliable  indicator  was  a  marked  decrease  in  relative  humidity,  even 

for  cloudless  conditions.  A  cloudless  example  is  given  in  Fig.  4  ,  which 

2 

shows  profiles  of  Cn,  temperature,  and  relative  humidity  for  1417  MST, 
9/7/84,  at  WSMR  for  which  a  z^  of  2100  m  was  determined.  A  temoerature 
profile  with  a  -l°C/100m  slope,  corresponding  to  adiabatic  conditions ( r) , 
is  drawn  near  the  measured  temperature  profile  for  comparison.  It  can 
be  noted  that  the  measured  profile  is  nearly  adiabatic  above  about  300  m, 
with  a  small  change  in  stability  above  2200  m.  Relative  humidity, 
however,  (1)  increases  gradually  from  19  to  26%  from  13  m  to  1000  m, 

(2)  increases  sharoly  to  a  maximum  of  47%  between  1000  and  2100  m  and 

(3)  decreases  sharply  to  about  20%  between  2100  and  2400  m.  A  secondary 

2 

maximum  of  Cp  can  also  be  noted  near  2000  m. 

Reasons  for  the  marked  changes  in  relative  humidity  (R.H.)  may  be 


(9 


26  53  40 

RELATIVE  HUMIDITY  «J 


47 


Fig.  4 


Profiles  of  C^,  temperature  and  relative  humidity 
for  1417  MS7,  9/7/84  at  White  Sands,  NM. 
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accounted  for  by  considering  its  definition  in  percent  as 

R.H.  =  ~  (100)  ,  (10) 

ws 

which  is  the  ratio  of  the  actual  mixing  ratio  (w)  to  the  mixing  ratio 
of  air  saturated  at  the  same  temperature  (w$).  The  actual  and  satura¬ 
tion  mixing  ratios  are  the  mass  of  water  vapor  per  unit  mass  of  dry  air 
and  saturated  air,  respectively.  The  value  of  w  depends  on  dew  point, 
which  is  a  single-valued  function  of  the  amount  of  water  vapor  per  unit 
volume,  and  w$  depends  on  temperature  only.  Dew  point  is  less  than 
temperature  except  in  completely  saturated  conditions,  such  as  dense 
fog  for  which  it  can  equal  temperature, and  RH  nl  100%.  Because  the 
temperature  change  was  small  at  the  height  corresponding  to  the 
change  in  relative  humidity,  ws  changed  very  little,  but  the  actual 
amount  of  water  vapor  (w)  must  have  decreased  sharply  above  2100  m  to 
cause  the  decrease  in  relative  humidity.  Physical  processes  responsi¬ 
ble  for  such  a  large  change  in  cloudless  conditions  involve  the  struc¬ 
ture  of  thermals  and  complex  thermodynamics  that  are  beyond  the  scope 
of  the  present  research. 

For  stable  stratification,  values  of  2n  were  detectable  as  the 
height  at  which  a  ground-based  temperature  increase  with  height  changed 
to  a  near-adiabatic  decrease  that  was  coincident  with  an  increase  in 
relative  humidity.  By  graphing  all  values  recorded  within  only  600  m 
above  ground  instead  of  average  values  up  to  3000  m,  it  was  possible 
to  obtain  adequate  resolution  to  allow  zn  to  be  reliably  determined 
for  most  soundings. 
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2.2  Rome  Air  Development  Center  (RADC)  data 

2 

Magnetic  tape  recordings  of  3 -minute  averages  of  (1)  Cn  measured 

-? /3 

with  an  optical  scintillometer  (m  ''  ),  (2)  air  temperature  (°C),  (3) 
dew  point  (°C),  (4)  wind  speed  (mph),  (5)  wind  direction  (°az.),  (6) 
station  pressure  (mb),  and  (7)  global  (direct  plus  diffuse  on  a  hori¬ 
zontal  surface)  solar  irradiance  (Langley  min"^)  we^e  processed  and 
analyzed  for  selected  days  in  1982.  The  measurements  were  made  by  Air 
Force  personnel  at  the  Advanced  Optical  Test  Facility  of  the  Rome  Air 
Development  Center  (RADC),  Verona,  NY  (Fragapane,  et  ai.»  1983). 

The  scintillometer  system  used  at  RADC  was  based  on  the  design  of 
Ochs  et  al.  (1977)  and  consisted  of  a  quartz-iodine  light  source  located 
118  m  from  a  receiver.  The  receiver  had  two  5-cm  diameter  apertures 

spaced  13.5  cm  apart.  Fluctuations  in  refractive  index  cause  variations 

2 

in  light  intensity  at  the  two  apertures  and  Cp  is  proportional  to  the 

2 

difference  in  light  intensity.  The  system  responded  to  values  of  Cn 
-17  -12  -2/3 

between  about  10  and  10  m  and  was  sensitive  primarily  to 

turbulent  refractive  index  scale  sizes  of  about  5  cm.  The  optical  path 

had  a  north-south  orientation  and,  like  the  meteorological  sensors,  was 

2 

at  a  height  of  2  m.  Voltage  signals  from  the  Cj^  and  meteorological 
sensors  were  sampled  50  times  each  minute,  digitized,  and  each  3-minute 
average  was  entered  onto  magnetic  tape. 

Hourly  weather  observations  for  1982  were  obtained  for  Griffiss 
Air  Force  Base,  located  23  km  from  the  test  facility,  and  were  used 
to  select  daytime  and  nighttime  periods  for  analysis.  Because  the 
largest  values  of  K  10"  m  '  )  occur  with  a  cloudless  sky  and  light 
wind,  the  criterion  for  selecting  periods  for  analysis  was  that  the 
observed  sky  condition  at  Griffiss  AFB  be  either  cloudless  or  that  it 
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have  only  scattered  clouds  (half  or  less  of  the  sky  covered  with  clouds) 
and  that  the  visibility  exceed  10  miles,  with  both  conditions  lasting 
for  at  least  2  hours.  Daytime  was  taken  to  be  from  one  hour  after  sun¬ 
rise  to  one  hour  before  sunset,  and  nighttime  was  taken  to  be  from  one 
hour  after  sunset  to  one  hour  before  sunrise. 

It  was  assumed  that  average  meteorological  conditions,  such  as 
surface  visibil ity,  humidity  and  the  type,  height  and  amount  of  clouds 
observed  at  Griffiss  APB  were  representative  of  those  at  the  test  site 
14  miles  away.  Limited  information  obtained  concerning  the  test  site, 
however,  indicated  that  it  was  located  at  a  lower  elevation  and  had  a 

greater  than  average  soil  moisture,  so  this  assumption  was  probably  not  ] 

1 

valid  for  certain  conditions.  For  several  cloudless  nights,  for  example, 
temperature,  dew  point  and  wind  conditions  were  conducive  to  the  forma¬ 
tion  of  ground  fog  at  Griffiss  APB  but  none  was  reported.  If  it  had 

occurred  at  the  test  site,  attenuation  of  the  scintillometer  light  beam 
2 

could  render  Cn  data  useless.  As  a  criterion  in  screening  data,  there¬ 
fore,  if  the  temperature-dew  point  difference  measured  at  the  test  site 
was  less  than  2°C,  the  data  were  not  used.  In  addition,  differences 
in  cloudiness  between  the  two  sites  could  also  occur  and  are  probably 

p 

responsible  for  some  of  the  scatter  in  the  relationships  of  Cn  to 
meteorological  variables.  The  chance  of  occurrence  of  different  cloudi¬ 
ness  conditions  at  the  two  locations  could  have  been  minimized  by  select- 
ing  periods  only  with  a  cloudless  sky  observed  at  Griffiss  AFB,  but  then 
(1)  the  sample  size  would  have  been  severely  reduced  and  (2)  important 
effects  of  variable  cloud  conditions  on  Cn,  such  as  the  large  changes 
that  can  occur  in  short  time  intervals  when  there  are  occasional  shadows 
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from  scattered  cumulus  clouds  on  snowless  ground,  would  not  have  been 
disclosed. 

Periods  from  winter  months  when  there  was  at  least  3"  of  snow 
and  the  air  temperature  was  above  0°C  were  also  selected  for  analysis, 

regardless  of  cloudiness.  As  discussed  below,  these  are  special  cases 

2 

that  can  produce  large  values  of  Cn  because  a  snow  surface  temperature 

remains  near  0°C,  regardless  of  how  much  the  air  temperature  exceeds  0°C. 

2 

C  for  daytime  conditions  over  snow  with  T  .  <  0°C,  on  the  other  hand, 

n  air 

/  -15  -2/3  1 

was  found  to  be  small  (£10  m  )  and  quite  steady.  As  a  result,  s 

! 

periods  with  these  conditions  were  not  selected  for  analysis.  j 

In  summary,  RADC  data  were  analyzed  separately  for  the  following  \ 

conditions: 

1)  Nighttime 

a)  No  snow  cover,  cloudless  or  scattered  clouds,  regardless 
of  air  temperature; 

b)  With  snow  cover,  cloudless  or  scattered  clouds,  but  res¬ 
tricted  to  Tair  <  0°C;  and 

c)  With  snow  cover,  regardless  of  cloudiness,  but  restricted 

to  T  •  >  0°C. 
air 

2)  Daytime 

a)  No  snow  cover,  cloudless  or  scattered  clouds,  regardless 
of  temperature;  and 

b)  With  snow  cover,  regardless  of  cloudiness  but  restricted 

to  T  •  >  0°C. 
air 

The  nighttime  periods  are  listed  in  Table  5  ,  together  with  average 

2 

cloud  amount  and  height  and  the  maximum  and  minimum  3-minute  values  of  C^, 
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TABLE  5  Maximum  and  minimum  values  of  Cn  and  relevant 
meteorological  variables  for  nighttime  periods 
analyzed  from  the  RADC  Optical  Test  Facility, 
for  a  2-meter  measurement  height 


CLOUD1  WIND  (mph)  TEMP  (°C)  DEW  PT  (°C) 


MO/ DAY  BEG.  END  MIN  MAX 


(AMT.)  HT.  MIN  MAX  MIN  MAX  MIN  MAX 


TABLE  5  (Continued.) 


1982  TIME  (GMT) 


MO/DAY  BEG.  END  MIN  MAX 


CLOUD1  WIND  (mph)  TEMP  (°C)  DEW  PT  (°C) 
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j  wind  speed,  temperature  and  dewpoint  for  each  period.  The  cloudiness 

i 

values  are  those  obtained  from  the  Griffiss  AFB  observations  and  are 

2 

representative  averages  for  each  period.  Values  of  less  than 
-17  -2/3 

1  x  10  m  '  were  not  used,  and  if  the  measured  wind  speed  was  less 
than  0.5  mph,  it  is  noted  with  a  single  asterisk.  Periods  are  listed 
for  43  nights  that  met  the  above  criteria  and  in  addition,  had  most 
equipment  operating.  After  final  screening,  4041  3-minute  values  com¬ 
prised  the  nighttime  RADC  data  set. 

The  daytime  periods  are  listed  in  Table  6  together  with  average 
cloud  amount  and  height,  estimates  of  the  top  of  the  mixed  layer  (z^) 
at  1800  GMT  with  the  Smith  (1977)  nomogram,  and  maximum  and  minimum 
values  of  Cn,  solar  irradiance  and  wind  speed.  Periods  from  20  days, 
or  about  213  hours  with  4267  3-minute  values,  comprise  the  daytime 
RADC  data  set. 


■J 

'i 
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o 

TABLE  6  Maximum  and  minimum  values  of  and  relevant 
meteorological  variables  for  daytime  periods 
analyzed  from  the  RADC  Optical  Test  Facility? 
for  a  2 -meter  measurement  height 

1982  TIME  (GMT)  C2x10"14(nf  2/3)  CLOUD1  SOLAR  (LY/MIN)  WIND  (mph)  Z^m) 

MO/DAY  BEG.  END  MIN  MAX  (AMT.)  HT.  MIN  MAX  MIN  MAX  18  GMT 


1115  2248  0.059 


1109  2251 


0.051 


1100  2257  0.052 

1042  2315  0.001 


1024  2333  0.054 


1024  2342 
1024  2345 


1027  2345 


0.004 


1027  2342  0.001 

1027  2342  0.002 

1024  2333  0.0002 


1057  2315  0.006 


9/11  1136  2221 

9/13  1136  2218 

10/4  1729  2139 

10/5  1203  2136 

11/17  1257  2036 


0.005 


1351  2024  0.004 


12/17  1533  2018 

12/18  1336  2024 


0.002 


0-4  40 

1- 5  40 
0-1  200 
0-5  250 
0-1  200 

2- 5  30 

1-4  30 

1-4  30 

3- 5  35 

0-4  40 

0-2  30 

1-4  35 

0-2  40 

0-3  40 


0-6  250 
0-6  100 
2-5  35 

1  200 


1-3  120 


*  14.9 

1.3  21.2 
2.9  12.7 

*  11.7 
5.2  16.5 
2.5  13.0 

1.4  14.9 

2.5  13.8 

*  10.1 


1.1  I  7.9 


7.3  10.3 


2.4  14.5 


*  <0.5  mph 

1  *  CLOUD  AMOUNT  IN  TENTHS  OF  SKY  COVERED;  CLOUD  HEIGHT  IN  HUNDREDS  OF  FEET 


mt 


2.3  NOAA/ERL  Boulder  Atmospheric  Observatory  (BAO)  data 

2  2 

Information  on  the  variation  of  Cn(Cy)  and  meteorological  variables 
between  10  m  and  300  m  above  ground  was  obtained  from  measurements  made 
at  the  Boulder  Atmospheric  Observatory  (BAO),  CO.  The  facility  is 
operated  by  the  NOAA/ERL  Wave  Propagation  Laboratory  and  is  located 
about  25  km  east  of  the  foothills  of  the  Rocky  Mountains  in  Colorado. 

The  profile  measurements  were  made  on  a  300-m  tower  there  and  are  unique 
in  their  type,  detail  and  format  for  analysis  (Kaimal  and  Gaynor,  1983). 

2 

The  BAO  tower  instrumentation  for  the  measurement  of  Cy  and  meteoro¬ 
logical  profiles  was  (1)  a  3-axis  sonic  anemometer  for  wind  component 
fluctuations,  (2)  a  propeller-vane  anemometer  for  average  wind  velocity, 
(3)  quartz  (slow  response)  and  platinum-wire  (fast  response)  temperature 
sensors  and  (4)  a  cooled  mirror  dew  point  hygrometer.  This  combination 
of  sensors  was  located  at  heights  of  10,  22,  50,  100,  150,  200,  250  and 
300  m  above  ground.  In  full  operation,  mean  and  fluctuation  measurements 
of  temperature,  wind  and  humidity  and  other  variables  were  automatically 
recorded  in  20-minute  segments.  Otherwise,  only  mean  quantities  were 
recorded.  Coincident  ground-based  measurements  included  solar  irra- 
diance  on  a  horizontal  surface  and  Cp  obtained  with  scintillometers  at 
a  height  of  4  m  in  three  450-m  optical  paths  in  the  shape  of  a  triangle. 
The  purpose  of  this  scintillometer  arrangement  was  to  enable  wind  con¬ 
vergence  characteristics  to  be  measured  optically  (Kjelaas  and  Ochs, 
1974). 

The  processing  and  analysis  of  BAO  data  were  carried  out  on  the 
University  of  Michigan  AMDAHL  470  V/8  computer.  The  tape-recorded  20- 
minute  values  of  mean  temperature,  dew  point,  wind  speed  and  temperature 
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fluctuation  spectral  data  for  the  8  heights  were  printed  out  and  graphed. 

2 

Values  of  Cy  for  each  height  were  calculated  from  spectral  densities  for 

_3 

22  frequencies  between  10  and  4  Hz.  With  graphs  of  successive  spectra, 

it  was  possible  to  determine  limits  of  the  inertial  subrange  f %om  the 

o 

-2/3  slope  relationship  (Kaimal,  1973).  Values  of  Cy  were  calculated 
with  the  spectral  density  corresponding  to  a  particular  frequency  in  the 
inertial  subrange  with  the  following  relationship: 

Cj  -  13.6  (n/lD2/3  n  ST(n)  ,  (n) 

o 

where  Sy(n)  is  the  spectral  density  of  temperature  fluctuations  (°C  Hz  ) 

at  n  =  2.03  Hz  and  U  is  the  corresponding  mean  20-minute  wind  speed 

2 

measured  at  each  height  of  Cy.  The  coefficient  13.6  results  *rom  the 

transformation  from  frequency  to  wave  number.  Profile?  of  fcn  C2  vs  tn  i 

-4 ;  3 

were  graphed  by  computer  and  compared  with  a  :  1  variation  for  corre- 

lation  with  meteorological  conditions.  Values  of  Cy  were  also  converted 
2 

to  Cn  with  eq.  (3)  for  comparison  with  corresponding  scintillometer 
2 

measurements  of  Cn-  Data  for  that  scintillometer  with  an  optical  path 
oriented  in  a  northeast-southwest  direction  were  used. 

Twenty-minute  data  for  daytime  and  nighttime  conditions  were  obtained 
for  (1)  11  days  from  a  month-long  intensive  field  experiment  conducted 
in  September  1978  called  Project  PHOENIX  (Hooke,  1979),  (2)  those 
periods  closest  to  the  launch  times  of  thermosondes  by  AFGL  personnel 
from  Table  Top  Mountain  near  Boulder  on  7  days  in  March  1983  and  (3) 

23  days  selected  from  June  through  September,  1983.  As  much  as  possible, 
the  data  were  selected  only  for  a  cloudless  sky  or  one  with  scattered 
clouds.  The  PHOENIX  periods  were  documented  with  detailed  weather 


i 


observations  and  other  meteorological  data,  but  for  the  remaining 
pjriods,  cloudiness  conditions  at  the  BAO  location  were  inferred  from 
weather  naps  and  hourly  observations  made  at  the  National  Weather  Ser¬ 
vice  Station  at  Denver's  Stapleton  Airport,  about  40  km  south  of  Boulder. 

For  most  days  conditions  were  similar,  but  it  is  likely  that  for  some 
periods,  the  sky  may  have  been  overcast  at  Denver  and  cloudless  at  BAO, 
or  vice  versa. 

The  daytime  periods  are  listed  in  Table  7  in  a  format  similar  to 
that  used  for  the  RADC  data.  There  were  227  sample  20-minute 
periods  for  the  1978  PHOENIX  daytime  data  and  325  samples  for  the  1983 
data.  The  Cn  values  are  scintillometer  measurements  and  the  values  of 
wind  speed  are  those  at  10  m  on  the  tower.  Also  shown  for  several  days 
in  the  September  1978  PHOENIX  data  are  values  of  2^  reported  by  Kaimal  (1982) 
from  measurements  with  various  active  and  passive  sensing  equipment. 

They  are  listed  with  values  estimated  from  the  Smith  (1977)  nomogram. 

For  the  1983  data,  only  values  estimated  with  the  Smith  nomogram  are 
available.  Nighttime  periods  and  relevant  data  for  the  BAO  data  set  are 
given  in  Table  8  .  A  total  of  360  20-minute  periods  comprise  the  night- 


Table  7  Values  of  Cn  and  relevant  meteorological 
variables  for  daytime  periods  analyzed 
from  the  Boulder  Atmospheric  Observatory 


TIME  (MST) 


BEG  END 


CLOUD2  SOLAR  (LY/MIN)  WINCr(ms"  )  Z.(m)(13MST) 


MIN  MAX  AMT  HT  MIN  MAX  MIN  MAX  MEA 


9*3 


0720 

0800 

1200 

0840 

0640 

0640 

0720 

1020 

1420 

0940 

0740 

0740 

0700 

0700 


a 


Table  7 


TIME  (MST) 


BEG  END 


AUG  1983 


31  0840  1300 
1340  1700 


SEP  1983 


1  I  0720  1140 
1220  1240 
1340  1420 
1520  1700 
5  I  0840  1700 


6  | 0640  1300 
1340  1700 

7  I  0720  1400 


1440  1700 
9  0820  1120 
1220  1700 


22  0640  0920 


1000  1100 
1140  1340 


Values  of  Cp  and  relevant  meteorological 
variables  for  daytime  periods  analyzed 
from  the  Boulder  Atmospheric  Observatory 
Tower  (continued) 


CLOUD^  SOLAR  (LY/MIN) 


MIN  MAX  AMT  HT  MIN  MAX  MIN  MAX  MEA 


0-4  8 


11.5  52.9  0 


23.3  91.9  0 

1.73  46.1  0 


1.14  30.2  0 


1.19  26.3  0 


Table  8 
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Values  of  and  relevant  meteorological 


variables  for  nighttime  periods  analyzed 
from  the  Boulder  Atmospheric  Observatory 
Tower 


TIME  (MST) 


AlO'-^m'^3)1 


CLOUD*1 


WIND3(ms"1 ) 


TEMP  (°C)  DEW  PT  (°C) 
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I 


2 

2.4  Compatibility  of  various  Cn  measurements 

2 

Because  estimates  of  C  were  made  from  measurements  with  various 

n 

optical  and  temperature  fluctuation  equipment  in  different  locations, 

2 

it  is  reasonable  to  expect  Cn  data  sets  that  are  internally  consistent 
but  somewhat  different  from  each  other  because  of  differences  in  (1} 
physical,  electrical  and  optical  characteristics  of  the  equipment  used 
and  (2)  thermal,  roughness  and  radiative  characteristics  of  the  measure¬ 
ment  sites  and  (3)  methods  of  data  acquisition  and  processing. 

As  pointed  out  by  Kunkel  and  Walters  (1982),  scintillometer  measure¬ 
ments  tend  to  converge  to  a  mean  very  rapidly  compared  with  point  measure¬ 
ments  made  with  either  a  single  probe  or  with  2  spaced  probes.  The  main 
reason  is  that  because  scintillometer  measurements  are  path  averaged,  a 
single  sample  in  time  is  already  an  ensemble  average.  In  addition,  the 
sensitivity  of  a  scintillometer  primarily  to  small  turbulent  eddy  sizes 
U  5  cm)  causes  the  statistical  uncertainty  to  decrease  rapidly  with  time 
because  effects  of  a  large  number  of  eddies  can  be  measured.  A  fast  con¬ 
vergence  of  these  measurements  allowed  Kunkel  and  Walters  to  measure 

rapid  changes  in  flow  characteristics. 

2 

With  regard  to  the  relationships  obtained  with  the  scintillometers 

at  RADC  and  BAO,  even  though  their  optical  and  electrical  characteristics 

were  very  similar,  their  different  physical  characteristics  alone  resulted 

2 

in  different  values  of  Cn  for  similar  meteorological  conditions.  The 
differences  could  have  been  caused  by  the  2-m  height  of  the  RADC  optical 
path  as  opposed  to  the  4-m  height  of  the  BAO  path  because  of  the  ex¬ 
ponential  decrease  of  with  height  in  the  first  few  meters.  For 


similar  meteorological  conditions  in  unstable  stratification,  for  example, 

2 

effects  on  Cn  of  the  2-m  path  height  at  RADC  and  the  4-m  height  at  BAO 

2  -4/3 

alone  should  cause  Cn  at  BAO  to  be  less  by  a  factor  between  about  (4/2) 

-2/3 

=  0.39  and  (4/2)  =  0.63,  depending  mainly  on  the  thermal  and  wind 

structure  in  the  first  few  meters.  Even  though  the  scintillometer  path 

length  is  118  m  at  RADC  and  400  m  at  BAO,  a  path  length  effect  alone 
2 

should  not  cause  Cn  at  BAO  to  be  significantly  greater  at  BAO  as  it  would 

be  if  measured  with  a  single  receiver  and  light  source  (Portman  et  al.,  1962). 

Quantifying  effects  of  these  inherent  physical  differences  separately 

was  not  possible  with  the  information  available,  especially  since  thermal, 

roughness,  and  radiative  characteristics  of  the  measurement  sites  were 

also  different.  An  indication  of  combined  effects  on  scintillometer 
2 

values  of  CR,  however,  can  be  noted  in  Fig.  5,  which  compares  data  from 
RADC  with  the  BAO  relationship  shown  in  Fig.  7.  It  shows  that  for  a  given 
solar  irradiance,  Cn  is  at  least  50%  greater  at  BAO  than  at  RADC.  Appa¬ 
rently,  effects  of  the  different  scintillometer  type  at  BAO  in  increasing  CR  sig 
nificantly  exceed  effects  of  the  higher  path  in  decreasing  it,  assuming 
that  average  effects  of  other  meteorological  variables  that  determine 

Cn  are  taken  into  account  by  the  large  sample  sizes. 

2 

Values  of  Cn  obtained  from  temperature  fluctuation  spectra  as  point 
measurements  can  also  be  expected  to  be  somewhat  different  from  the 
integral -path  measurements  with  scintillometers,  depending  on  meteorolo¬ 
gical  conditions.  A  conflicting  result  that  occurred  in  the  BAO  data  is 
that  for  a  given  solar  irradiance,  CR  obtained  from  temperature  spectra 
at  10  m  was  about  38%  greater  than  Cn  from  the  scintillometer  measure¬ 
ments  at  4  m,  whereas  it  would  be  expected  to  be  less  by  a  factor  between 
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about  (10/4)“4//3  =  0.29  and  (10/4 )~2//3  =  0.54.  Reasons  for  this  differ¬ 
ence  were  not  determined,  but  in  addition  to  those  given  above,  they 
involve  questions  regarding  the  accuracy  of  spectra  at  low  wind  speeds, 
assumptions  regarding  the  isotropy  of  turbulence  and  possible  effects 
of  a  tower  structure  on  temperature  spectra  for  certain  wind  directions. 
For  these  reasons,  values  of  obtained  from  BAO  temperature  spectra 
were  analyzed  independently  of  those  obtained  from  the  scintillometer. 

A  comparison  was  also  made  between  scintillometer  and  temperature 
probe  estimates  of  Cn  at  heights  of  4  and  14  m  for  2  cloudless  days  (7 
and  8  Sept.  1984)  at  White  Sands,  NM.  The  4-m  and  14-m  scintillometers 

had  path  lengths  of  400  m  and  1  km,  respectively.  With  the  temperature 

? 

difference  sensors,  the  temperature  structure  parameter  C-j.  was  determined 
with  the  expression: 

C2  =  UT}2  R‘2/3  ,  (12) 

where  aT  is  the  rms  value  of  the  temperature  difference  and  R  is  the 
distance  between  sensors,  assuming  isotropy  of  turbulence  and  that  R  is 
within  the  inertial  subrange.  A  listing  of  daytime  values  of  meteorolo¬ 
gical  variables  and  for  every  hour  on  the  hour  for  the  2  days  is 
given  in  Table  9.  Occasional  large  differences  are  evident  for  a  given 
height  that  could  be  due  to  the  factors  discussed  above  in  addition  to 
the  visual-manual  method  used  to  abstract  the  values,  but  in  general, 
acceptable  agreement  for  these  and  other  days  was  found. 


Table  9  Values  of  meteorological  variables  and  Cn 
on  the  hour  for  2  cloudless  days  at  White 
Sands,  NM 


DATE 

TIKE 

(MST) 

T  ~T  ^ 
SFC  ‘AIR 

°C 

DIR. 

WIND2 
SP(ms“^ ) 

SOLAR3 

Wnf2 

c>' 
SC  I  NT.4 

4m  14m 

■’V2/3) 

at  probes5 

4m  14m 

9/7/84 

07 

3.3 

360 

0.7 

200 

5.7 

1 

14 

40 

08 

8.9 

280 

0.7 

413 

15 

0.4 

9.7 

7 

09 

11.1 

290 

2.8 

640 

19 

5 

7 

3.1 

10 

13.9 

360 

2.2 

813 

29 

6.7 

11 

6.; 

11 

18.9 

330 

1.8 

884 

38 

7.4 

35 

9.( 

12 

19.4 

310 

3.8 

917 

38 

11 

37 

8.' 

13 

18.3 

320 

4.5 

891 

35 

11 

34 

10 

14 

17.2 

300 

3.5 

815 

38 

8.5 

30 

7 

15 

15.0 

310 

3.6 

660 

28 

6.7 

24 

6.' 

16 

7.8 

300 

2.6 

451 

17 

3.4 

12 

3.1 

17 

1.7 

290 

2.3 

227 

4.8 

0.8 

1.7 

O.l 

9/8/84 

| 

07 

3.3 

360 

0.6 

180 

5.2 

1.2 

8 

6.: 

08 

8.3 

360 

0.4 

412 

14 

2 

14 

13 

09 

12.8 

360 

1.4 

635 

25 

5.3 

18 

17 

10 

16.7 

350 

1.4 

804 

35 

6.3 

42 

42 

11 

20.0 

230 

1.3 

892 

45 

7.5 

41 

12 

12 

21.1 

210 

1.3 

917 

50 

10 

45 

10 

13 

20.6 

220 

1.6 

895 

48 

8.6 

45 

17 

14 

19.4 

360 

1.6 

809 

45 

7.2 

39 

15 

15 

15.6 

300 

2.2 

667 

28 

6.7 

28 

11 

16 

7.8 

290 

1 .7 

457 

15 

3.3 

13 

3. 

17 

1.7 

310 

2.5 

224 

2.4 

0.7 

1 

0. 

Tair  is  average  of  temperatures  at  0.5,  1,  2,  4  m  estimated  from  curve  drawn 
through  graphed  diurnal  variation,  tower  2(HIDL). 

2 

Hind  velocity  at  4m, tower  1  (HIDL). 

3 

Global  solar  irradiance  at  14  m,  tower  1  (HIDL). 

4 

Scintillometer  path  length  is  400  m  at  4  rn  and  1  km  at  14  m  (HIDL). 

^  Horizontal  temperature  difference  probes,  tower  1  (HIDL). 


-  48  - 


3,  Data  analysis 

The  approach  in  the  analysis  was  to  develop  and  test  relationships 

2  2 

that  will  (1)  provide  estimates  of  C^(Cy)  for  one  height  in  the  atmos¬ 
phere's  first  few  meters,  (2)  estimate  the  depth  of  the  daytime  con¬ 
vective  boundary  layer  with  solar  radiation,  cloud  and  wind  speed 
information  and  (3)  apply  an  expression  similar  to  that  of  Kukharets 

p 

and  Tsvang  (1980)  to  obtain  Cn(z).  It  is  assumed  that  because  the 
2 

magnitude  of  Cn  near  the  ground  is  determined  by  thermal  and  mechanical 

2 

turbulence  processes  linked  by  similarity  principles,  the  value  of  Cn 
at  one  height  can,  in  most  cases,  be  a  reliable  indication  of  its  pro¬ 
file  above  or  below  that  height.  As  an  extreme  example,  with  adiabatic 

2 

conditions  throughout  the  boundary  layer,  minimum  values  of  can  be 

2 

expected.  Not  until  Cn  reaches  a  certain  value  at  a  given  height 
should  it  reflect  the  thermal  and  wind  structure  described  by  similarity 
theory  that  leads  to  a  particular  variation  with  height. 

3.1  cn(Cj)  in  unstable  stratification 

Determining  effects  of  solar  radiation  was  emphasized  for  the 
following  reasons:  (1)  the  amount  of  direct  and  diffuse  solar  radia¬ 
tion  reaching  the  ground  is  the  primary  determinant  of  how  much  heating 
of  the  ground  takes  place  and  how  unstable  the  lapse  rate  in  the  surface 
layer  will  become  as  it  is  acted  on  by  wind  to  produce  Cn(Cy),  (2)  it 
can  be  estimated  for  time  intervals  as  short  as  an  hour  if  location 
(latitude,  longitude),  date,  time  of  day,  and  cloud  type,  height  and 
amount  are  known  (Kasten  and  Czeplak,  1980)  and  (Turner  and  Mujahid, 
1984),  and  (3)  it  is  implicitly  important  in  determinations  of  z^,  as 
shown  in  Fig.  2  (Smith,  1977). 
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In  this  section,  relationships  among  these  variables  are  deve¬ 
loped  for  the  RADC  and  BAO  data  for  a  single  height  and  tested  with 

2 

hourly  values  of  abstracted  from  graphs  for  17  days  for  CLEAR  I. 

2 

Also,  measured  boundary  layer  profiles  of  Cp  for  CLEAR  I  and  II  are 

2 

compared  wich  profiles  calculated  with  the  K-T  model,  using  Cn  measured 
at  a  reference  height  and  values  of  obtained  from  temperature  and 
humidity  profiles. 

3.1.1  Single-height  relationships 

For  the  RADC  data  3-minute  values  of  Cp  were  graphed  against 
solar  irradiance  for  individual  days  shown  in  Table  6.  For  most  days, 

Cn  increased  linearly  with  solar  irradiance  from  an  early  morning 
minimum  just  after  sunrise  to  a  midday  maximum  near  solar  noon  and  then 
decreased  linearly  to  a  late  afternoon  minimum.  The  afternoon  minima 
of  Cn  occurred  at  higher  values  of  solar  irradiance  than  the  morning 
minima  and  are  explainable  in  terms  of  energy  exchange  processes  and 
the  time  of  occurrence  of  a  near  adiabatic  temperature  profile  after 
sunrise  and  before  sunset. 

For  most  of  the  RADC  days  analyzed,  it  was  found  that  for  a 

given  solar  irradiance  either  side  of  solar  noon,  mid-afternoon  values 

of  Cn  were  commonly  half  those  in  mid-morning.  An  example  is  shown  in 

Fig.  6,  which  is  a  computer  graph  of  values  for  15  May  1982,  a  day  with 

about  2/10  cloudiness  and  a  wind  speed  of  3  ms"^ .  It  can  be  noted  that 

a  morning  increase  of  Cn  to  a  maximum  has  values  that  are  about  twice 

those  of  an  afternoon  decrease  to  a  minimum.  Although  less  frequent 

and  less  pronounced,  this  difference  also  occurred  for  a  few  days  in 

2 

the  BAO  scintillometer  data  and  was  evident  in  a  sample  C  recording  for 


1 


.31  .46  .62  .77  .93  1.08  1.24 

GLOBAL  SOLAR  IRRADIANCE  (LY  MIN1) 


Fig.  6  Cn  versus  global  solar  irradiance  for  15  Hay  1982 
for  the  Rome  Air  Development  Center,  NY  (RADC). 
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one  day  for  Table  Mountain,  CO  (Ochs,  personal  cummunication) .  It 
was  not  evident,  however,  in  hourly  values  abstracted  for  7  cloudless 
days  in  the  CLEAR  I  data  for  White  Sands,  NM. 

If  the  morning  to  afternoon  differences  at  RADC  are  accepted 
as  being  due  to  natural  causes  and  not  to  the  equipment  used,  a  rea¬ 
sonable  explanation  is  possible  that  involves  the  availability  of  soil 
moisture  for  evaporation.  If  adequate  soil  moisture  is  available, 
differences  in  energy  exchange  processes  involving  latent  and  sensible 
heat  transfer,  net  radiation  and  soil  heat  transfer  from  morning  to 
afternoon  result  in  greater  vertical  temperature  differences  before 
solar  noon  than  afterward,  "external"  conditions  remaining  steady. 

Such  a  soil  moisture  condition  could  be  expected  most  of  the  time  in 
the  moist,  grass-covered  RADC  test  site,  occasionally  at  BAO  following 
rain,  but  seldom  at  White  Sands,  where  water  content  (%  by  weight)  in 
September,  1981,  was  found  to  range  from  near  zero  to  about  7 %  (Smith, 
1984).  Profile  and  other  measurements  are  necessary  to  test  the  expla¬ 
nation  given  above. 


Because  (1)  morning  to  afternoon  differences  in  Cn  were  rare 


and  not  as  pronounced  at  BAO  and  (2)  Cn  was  measured  at  4  m  and  the 


wind  speed  was  measured  at  a  standard  height  of  10  m,  the  BAO  data  were 


used  to  develop  the  final  regression  equation  involving  Cn,  solar 


irradiance  and  wind  speed.  Fig.  7  shows  average  values  of  Cn  versus 


solar  irradiance  (RAD)  for  the  BAO  data  that  were  used  to  develop  the 
regression  equation.  Standard  deviations  are  shown  as  vertical  line 
segements  and  sample  sizes  are  given  below  each  average. 


In  Fig.  7  a  nearly  linear  increase  of  Cf1  from  small  values 


at  irradiances  less  than  0.3  ly  min  to  maximum  values  at  irradiances 
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Fig.  7  Cn  versus  global  solar  irradiance  for  the  BAO 
scintillometer  data. 


i 
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near  1.25  ly  min  can  be  noted,  with  a  decrease  in  Cn  for  irradiances 

between  about  1.25  and  1.47  ly  min"1.  For  RAD  £  0.3  ly  min"1,  the 

small  values  of  Cn  are  representative  of  the  morning  and  afternoon 

transitions  in  aT /az,  as  Cn  approaches  the  noise  level  of  a  scintillo- 

<  <  -1 

meter  system.  For  0.3  -v  RAD  ^  1.25  ly  min  ,  a  linear  increase  in 
Cn  reflects  the  increase  in  aT/az  with  increased  solar  heating  of  the 
surface.  For  1.25  £  RAD  £  1.5  ly  min"1  a  decrease  in  Cn  is  more 
difficult  to  interpret  without  measurements  of  aT/az  because  the  reason 
for  Cn  leveling  off  and  then  decreasing  could  be  meteorological,  within 
the  Cn  measurement  system  or  a  combination  of  both.  The  last  2  data 
points  with  fewer  samples,  for  example,  are  for  very  high  irradiance 
values  that  are  likely  comprised  of  brief  increases  that  are  commonly 
caused  by  multiple  reflections  from  well -developed  cumulus  clouds. 

In  addition,  with  the  large  values  of  aT/az  that  probably  occur  in  this 
region,  light  source-to-receiver  alignment  problems  could  be  caused  by 
large  mean  refraction  in  the  optical  path  (see  also  RADC  data,  Fig.  5). 

Wind  speed  effects  were  determined  from  graphs  of  Cn  vs  10-m 
wind  speed  (u^ q )  for  various  categories  of  solar  irradiance,  from 
which  it  was  found  that  for  irradiances  greater  than  about  0.5  ly  min'1, 
Cn  increased  with  u^g,  reached  a  maximum  near  u^g  3  2-3  ms"1  and  gra¬ 
dually  decreased  at  higher  speeds. 

The  regression  equation  involving  C  ,  solar  irradiance  (RAD) 

2  3 

and  wind  speed  (WS)  was  developed  *>n  terms  of  RAD,  RAD  ,  and  RAD  and 
2  3 

WS,  US  ,  and  WS  .  Third-order  polynomials  were  chosen  for  these 
variables  to  account  for  (1)  the  two  inflection  points  for  Cn  vs  RAD 
shown  in  Fig.  7  and  (2)  the  one  inflection  point  for  the  C  vs  WS 
relationship.  The  equation  for  the  8A0  data  is: 
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Cp  =  .12724  x  10"6  +  .17314  x  10~6(RAD)  -  .33899  x  10'7(WS) 

+  .62238  x  10"6(RAD)2  +  .10264  x  10‘7(WS)2  (13) 


-  .39727  x  10"6(RAD}3  -  .53824  x  10“9(WS)3 


where  RAD  is  solar  irradiance  in  ly  min"^,  WS  wind  speed  in  ms"\  and 


I 


C  in  m 
n 


•1/3 


The  equation  is  valid  for  irradiances  between  0.07  and  1.49 

ly  mirf^  and  for  10-m  wind  speed  from  about  0.3  to  19  ms"^ .  It  can 

be  noted  that  a  value  of  Cn  <  0  is  possible  for  certain  values  of 

irradiance  and  wind  speed.  This  is  a  consequence  of  the  method  of 

-8  -1/3 

regression.  A  value  of  Cp  on  the  order  of  10“  m"  can  be  assumed 

for  a  calculated  C  <  0. 

n 

3.1.2  Boundary  layer  profiles 

2 

Fig.  8  shows  average  profiles  of  Cy  for  8  heights  between  10  m 

and  300  m  for  5  daytime  periods  of  the  1978  PHOENIX  experiment.  Because 

heignt  is  on  the  ordinate,  the  line  with  the  -3/4  slope  shown  for  refer- 

2  -4/  3 

ence  actually  corresponds  to  Cy  «  z  .  It  can  be  noted  that  except 
2 

for  values  of  Cy  at.  heights  between  200  and  300  m  on  26  September,  a 
-4/t 

z  relationship  is  followed  quite  closely.  For  all  the  profiles  in 

Fig.  8,  measurements  of  z  •  were  made  with  various  methods  (Kaimal  et  al., 

» 

1980  and  1982)  and  showed  that  z.  was  above  the  300-m  height.  On  26  Sep¬ 
tember,  however,  it  was  not  wel 1 -defined  and  could  not  be  detected  as 
an  increase  in  temperature  with  conventional  radiosonde  temperature 
measurements,  z.,  as  a  capping  inversion  and  having  other  properties 
detectable  by  all  methods,  was  most  pronounced  on  21  September,  when 

it  increased  from  500  m  to  approximately  850  m  during  the  period  repre- 
2 

sented  by  the  Cy  profile. 


\ 


Average  profiles  of  Cj  for  5  daytime  periods  from  the 
NOAA  Boulder  Atmospheric  Observatory  for  the  1978 
PHOENIX  exoerTment. 
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A  comparison  was  made  between  available  measurements  of 
reported  by  Kaimal  and  estimates  of  z.  from  the  Smith  (1977)  nomogram 
shown  in  Fig.  2.  Results  for  21,  22,  26  and  27  September  are  shown  in 
Table  10. 

Table  10  Measured  and  estimated  values  of  Zi  for  1300  MST  for 
BAO  data  from  PROJECT  PHOENIX 


Date 

>t.  1978) 

Cloudiness 

u10 

(ms-1) 

zi 
Meas . 

(m) 

Calc. 

21 

clear 

1.0 

850 

750 

22 

2/10  cirrus 

1.3 

480 

650 

26 

clear 

3.0 

650 

900 

27 

clear 

1.5 

680 

800 

2 

Comparisons  were  made  between  proxies  of  Cn  obtained  with 

2 

the  K-T  model  and  values  of  Cn  measured  with  thermosondes  during  CLEAR  I 
&  II  and  height-averaged  with  the  method  given  in  Table  2.  The  results 
of  the  comparisons  made  for  soundings  with  convective  conditions  and 
less  than  0.5  cloudiness  are  shown  in  Fig.  9a  through  9i.  In  each 
figure,  the  coordinates  are  logarithmic  and  the  abscissa  scale  for 
log  C“  is  the  same  ( ^  10  to  10  ).  The  numbers  or  the  ordinate  are 

actual  heights  that  correspond  to  the  logarithmic  height  scaling  by  com- 

9 

puter.  The  solid  line  is  the  Cn  profile  calculated  with  the  K-T  model 
and  the  dots  are  measured  values.  Values  of  CJ^(z)  were  calculated  with 
the  K-T  model  equation  as  modified  by  Murphy,  Dewan,  and  Sheldon  (1985), 
hereafter  referred  to  as  M0S,  to  be  applicable  to  any  reference  height, 
;r,  instead  of  :<r  =  0.1  as  in  Eq.  9:  The  expression  for  C^(z)  used  in 
the  calculations  was : 


i 

l 


1348  MST 
9/4/84 


t-  296 

o 

x 


x 

o 

3  94 


-17.1  -15.9 

LOG  C* 


.  J*  *  ^ 

•  • 


\  f  * 


i  319 

l  183 


1319  MST 
9/5/34 


LOG  c: 


Fig.  9a  Comparison  of  Cfj  profile  calculated  with  the  Kukharets  and 
Tsvang  model  ana  measured  values  for  unstable  stratification 
at  White  Sands,  NM  for  1348  MST,  9/4/84  and  1319  MST,  9/5/84 

Measured  values  are  height-averaged  as  m  Table  2. 


LOG  HGT,  HGT (m)  M  LOG  HGT;  HGT  (m) 


.  (14) 


cfc) 


C^0> 


_4/  -12(z/zrl.l)‘ 

0.046  (z/z.)  +  0.6  exp 

0.046  (z0/zi)-4/3 


where  Zg  is  a  reference  height,  taken  here  to  be  the  height  of  lowest 

measurement  of  Cn,  T  and  R.H.  for  each  sounding  (excluding  occasional 

"surface"  measurements),  z  is  height,  and  z..  is  the  top  of  the  mixed 

layer  obtained  from  profiles  of  temperature  and  relative  humidity. 

Actual  values  of  z^  and  Zg  for  each  sounding  are  given  in  Table  4  . 

Reasonably  good  agreement  between  the  K-T  model  results  and 
2 

measurements  of  Cn  is  evident,  including  the  region  with  the  secondary 
2  2 

maximum  of  C  near  z..  The  use  of  the  C  value  at  the  lowest  height  of 
n  i  n 

measurement  as  a  reference  value  for  the  model  profiles  worked  well  ex¬ 
cept  for  the  profile  on  1441  MST,  3/5/85.  For  this  profile,  the  re- 
9%  -14  -2/3 

ference  value  of  Cn  ^  4x10  m  for  Zg  =  44  m  is  unusually  large 

for  44  m  compared  to  values  at  the  same  height  for  other  soundings  and 

caused  all  calculated  values  to  be  larger  than  actual.  There  was  no 

2 

clear  evidence  that  the  distribution  of  C  in  the  region  of  z,  was 

n  ^ 

narrower  than  the  K-T  model  prediction  as  found  by  MDS,  at  least  not 
enough  narrower  to  support  their  value  of  k3  =  580  instead  of  12.  Some¬ 
what  different  results,  however,  could  be  expected  because  (1)  differ¬ 
ent  averaging  techniques  were  applied  by  MDS  to  the  C3  data,  and  (2)  in 
the  work  here,  the  values  of  z^  that  were  used  in  the  calculations  were 
also  probably  somewhat  different. 

The  results  support  other  findings  that  a  knowledge  of  at 
zQ  is  sufficient  information,  in  most  cases,  to  enable  a  profile  of 

in  a  C8L  to  be  described.  With  the  single-height  relationships  involving 
2 

Cn,  solar  irradiance  and  wind  speed  discussed  in  the  previous  section, 

2 

estimates  of  Cn  at  Zg  can  be  made  that  lead  to  boundary  layer  profiles 
with  the  MDS  modification  of  the  K-T  model. 


2  2 

3.2  Cn(Cy)  in  stable  strat ification 

3.2.1  Single-height  relationships 

In  contrast  to  the  combined  effects  of  positive  buoyancy 

2  2 

and  '  jrtical  wind  shear  that  enhance  turbulence  and  C  (CT)  in  a  convec- 

n  T 

tive  boundary  layer,  negative  buoyancy  and  viscosity  in  a  thermally 
stable  boundary  (SBL)  act  to  suppress  turbulence.  In  a  SBL,  mechanical 
turbulence  is  maintained  by  vertical  wind  shear  alone,  which  acts  on 
an  increase  in  temperature  with  height  to  produce  Cn(Cy).  At  the  air- 
ground  interface,  there  is  usually  a  nocturnal  loss  of  heat  from  the 
surface  by  infrared  radiation  and  a  gain  of  heat  by  (1)  a  downward 
transfer  of  infrared  radiation  and  sensible  and  latent  heat  from  air 
to  ground  and  (2)  an  upward  transfer  of  heat  by  conduction  through  the 
soil.  If  the  net  heat  exchange  is  negative,  the  surface  becomes  colder 
than  the  air  and  temperature  increases  with  height. 

Because  of  their  pronounced  effects  on  radiative  and  turbu¬ 
lent  exchange  processes  in  the  boundary  layer  regardless  of  time  of  day, 
cloudiness  and  wind  speed  are  the  primary  standard  meteorological  varia¬ 
bles  that  determine  nocturnal  temperature  and  wind  structure,  and,  there- 
2  2 

fore,  values  of  Cn(Cy).  As  clouds  become  fewer  and/or  higher  (colder),  for 
example,  the  net  radiative  exchange  and  surface  temperature  decrease  and 
air  tanperature  increases  with  height.  For  a  completely  cloudless  sky, 
radiative  loss  is  optimum  and  wind  speed  (shear),  through  its  mixing 
action,  becomes  the  controlling  variable  of  the  magnitude  of  the  inver¬ 
sion  in  the  first  few  meters.  In  the  absence  of  large-scale  advsction 
of  air  that  is  much  colder  than  the  surface,  the  higher  the  wind  speed, 
the  smaller  the  inversion.  For  nocturnal  conditions  with  radiative  loss. 
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2  2 

therefore,  because  Cp(Cy)  is  caused  by  wind  (shear)  acting  on  a  vertical 

2  2 

temperature  gradient,  it  was  reasonable  to  expect  that  values  of  Cn(Cy) 
at  one  height  as  well  as  their  profiles  would  be  relatable  to  wind 
speed.  By  confining  the  analysis  to  periods  with  a  sky  condition  that 
was  cloudless  or  scattered  ( <_  0. 5  cloud-covered),  furthermore,  a  large 
radiative  loss  leading  to  inversion  formation  could  be  expected,  leaving 
wind  speed  alone  as  the  most  likely  describer  of  Cn(Cy)  characteristics. 

In  line  with  this  reasoning,  the  4041  3-minute  values  com¬ 
prising  the  RADC  nighttime  set  of  data  for  the  periods  listed  in  Table  5 

were  grouped  into  1  mph  categories  of  the  wind  speed  at  2m,  and  the  average 
2 

values  of  Cn  and  wind  speed  were  obtained  in  each  category.  Categories 
with  fewer  than  30  samples  were  not  retained,  which  meant  that  there 
were  no  wind  speed  categories  with  sufficient  samples  greater  than  7  mph. 
The  results  of  the  grouping  are  shown  in  linear  coordinates  in  Fig.  10. 


Values  of  x  10'^m"2,/3 


are  on  the  ordinate  and  wind  speed  in  mph  is 


on  the  abscissa.  It  can  be  noted  that  CT  increases  from  a  minimum  of 
-14  -2/3 

0.5x10  m  at  a  wind  speed  of  *v  0.3  mph  to  a  maximum  of  about  6.2  x 

1Q*1V2/3  at  ^  4.3  mph  and  decreases  sharply  at  higher  speeds.  Limited 
2 

Cn  data  obtained  at  higher  speeds  show  that  values  become  progressively 

-17  -2/3 

smaller  as  wind  speed  increases,  approaching  about  10  in  '  . 

The  data  shown  in  Fig. 10  were  fitted  with  a  normal  curve 

shown  by  the  solid  line,  It  describes  the  RADC  relationship  quite  well 

for  nighttime  conditions  in  the  fora:  0 

(u?-4.3r 

cj  =  6.2  x  10'14  exp  2(1>8)2  ,  (15) 

1 4  2 

where  6.2x10*  H  is  the  peak  value  of  C^  occurring  at  a  2m  wind  speed  (u2) 


m 


1.1  2.2  3.3  4.4  5.5  6.6  7.7  8.8  10.0 

2-m  WIND  SPEED  (mph) 


Fig.  10 


2 

C  versus  2-n  wind  speed  on  linear  scales  for  nighttime 
conditions  at  RADC.  The  numbers  are  sample  sizes. 
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of  4.3  mph,  and  1.8  mph  is  the  standard  deviation  of  for  the  normal  curve; 
about  68%  of  the  area  under  the  curve  is  between  wind  speeds  of  2.5  and 

6.1  mph. 

Another  portrayal  of  the  RADC  data  is  given  in  Fig. 11  to 

illustrate  sample  sizes  and  standard  deviations  in  the  relationship.  It 

o 

has  on  a  logarithmic  ordinate  scale  and  wind  speed  on  a  linear  abscissa 
scale.  Sample  sizes  are  given  alongside  each  mean  and  the  vertical  lines 
are  standard  deviations.  It  can  be  noted  that  sample  sizes  range  from 
1970  for  u 2  =  3  mph  to  62  for  =  6.4  mph  and  that  the  standard  devia¬ 
tions  are  much  larger  for  u 2  <  3  mph,  indicating  turbulence  intermittency. 

Similar  results  were  obtained  with  the  scintillometer  data 
for  the  BAO  facility.  In  the  BAO  data,  however,  optical  Cp  was  measured 
at  a  height  of  4  m,  but  the  lowest  height  of  wind  speed  measurement  on 
the  300-m  tower  was  10  m.  Averages  of  Cn  (4  m)  and  u^Q  in  successive 
1  ms-1  categories  of  u^Q  are  shown  in  Fig. 12  together  with  sample  sizes 

and  a  normal  curve  estimated  to  fit  C  (4  m)  vs  u-jq.  It  can  be  noted 

o 

that  compared  to  the  RADC  results,  which  show  a  peak  of  (2  m)  = 

6.2  x  10-14m-2/3  at  =  4.3  mph  (1,9  ms’1),  the  BAO  data  have  a  peak  of 

&  (4  m)  =  12.8  x  lO^nf^3  near  u^Q  =  4.0  ms’1  (%  9  mph).  The  BAO  relation¬ 
ship,  fitted  with  the  equation  given  below,  is  the  nighttime  single-height 
o 

model  for  used  in  the  profile  calculations  in  Section  4.2. 

(u10-4.0)2 

cjj  =  12.8  x  10'Uexp2(2.3)2  ,  (16) 

-14  —2 

where  12,8  x  10  is  the  peak  value  of  (4  m)  occurring  at  a  10  m 

wind  speed  of  4.0  ms’1,  and  2.3  ms’1  is  the  standard  deviation  of  u1Q. 

2 

The  behavior  of  C  shown  in  Fig, 10  and  Fig. 12  is  caused  by 
n 

the  interaction  of  wind  speed  (shear)  with  the  thermally  stable  vertical 
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temperature  gradient.  With  few  or  no  clouds  and  no  detectable  wind  speed, 
the  largest  inversions  are  possible.  Because  vertical  motions  are  sup¬ 
pressed  by  the  large  negative  buoyancy  and  viscosity  associated  with  such 
inversions,  air  motion  is  mostly  wave-like  and  laminar,  with  slow  horizon¬ 
tal  meandering.  Visible  evidence  of  such  motion  can  be  seen  in  the  beha¬ 
vior  of  smoke  plumes  on  a  cloudless  and  windless  night.  For  such  a  con¬ 
dition,  observations  show  that  although  mean  optical  refraction  is  very 

2 

large,  optical  turbulence  itself  can  be  very  small,  with  Cn  values  occa¬ 
sionally  approaching  those  expected  for  adiabatic  conditions.  Such  small 
values,  however,  are  usually  only  temporary.  Abrupt  increases  occur  in¬ 
termittently  that  are  due  to 'bursts"  of  turbulence.  In  the  RADC  data, for 
example,  on  some  cloudless  nights  with  zero  wind  recorded  because  it  was 
less  than  an  anemometers  starting  speed,  the  3-minute  values  of  occa¬ 
sionally  increased  2  to  3  decades  in  short  (-v  10  minute)  periods  and 
lasted  for  about  10  minutes  before  decreasing.  Pronounced  effects  on 
optical  scintillation  and  visual  resolution  for  this  condition  occur 
that  have  been  measured  and  described  by  Portman  et  al .  ( 1 962)  and  by 
Ryznar  (1963).  Dewan  (1982)  presents  and  discusses  reasons  for  such 
behavior  in  terms  of  a  wave  interaction  and  cascade  processes. 

As  RADC  wind  speed  at  2  m  increased  from  nearly  calm  to  about 

4  mph,  the  magnitude  of  the  inversion  decreased  somewhat,  but  more  im- 

2 

portaritly,  vertical  motions  increased,  causing  a  large  increase  in  Cn. 

The  increase  was  more  evident  in  the  RADC  data  than  in  the  BAO  data 

because  3-minute  values,  rather  than  20-minute  values  were  available. 

If  conditions  were  conducive  to  a  steady  wind  speed  of  about  3-4  mph  at  RADC, 
2 

Cn  also  remained  steady  at  large  values. 
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As  wind  speed  continued  to  increase,  vertical  motions  also 
continued  to  increase,  but  the  inversion  magnitude  decreased.  As  shown 

O 

in  Fiq.10  and  12,  C  also  decreased.  The  decrease  in  inversion  magnitude  is 
due  primarily  to  enhanced  mixing  and  downward  turbulent  heat  transfer 
that  inhibit  cooling  of  the  surface  by  radiative  loss.  A  wind  speed 
was  eventually  reached  when  (1)  the  heat  lost  by  the  surface 

was  insufficient  to  keep  it  colder  than  the  air,  (2)  the  vertical 
temperature  gradient  probably  approached  adiabatic,  and  (3)  Cp 
decreased  to  steady  minimum  values  that  were  near  the  measurement  system's 
noise  level,  <v  3  ms"^  for  RADC  and  7  ms  ^  for  BAO. 

It  should  be  pointed  out  that  in  the  absence  of  marked  ad- 
vection  of  air  with  different  properties,  the  wind  speed  at  which  Cn 
minima  are  reached  are  greater  for  a  daytime  cloudless  condition  than 
for  a  nighttime  one.  The  daytime  gain  of  heat  by  solar  radiation 
at  the  surface  much  exceeds  that  lost  by  turbulent  and  radiative  heat 
transfer  and  keeps  the  surface  significantly  warmer  than  the  air  at 

much  higher  wind  speeds.  As  a  result,  values  of  Cn  with  2-m  wind 

-1  -13  -2/3 

speeds  between  about  3  and  5  ms  ,  for  example,  were  near  10  m 

_  I  r  _  p  f  n 

on  a  sunny  mid-day,  but  approached  10  m'  on  a  cloudless  night. 

Wintertime  conditions  with  a  snow  cover  and  with  air 

2 

temperatures  above  0°C  can  cause  large  values  of  Cn  in  stable  strati¬ 
fication  that  are  not  describable  in  terms  of  those  cloudiness,  solar 
irradiance  and  wind  speed  relationships  that  apply  to  conditions  with 

Tair  <  fr°r  Tair  <  ^  over  unifonnly  fresh  deep  snow  (>  10  cm) 

and  few  or  no  clouds,  the  average  temperature  of  a  snow  surface  (1) 
is  usually  colder  than  the  air  at  night  due  to  infrared  radiation  loss. 
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(2)  is  nearly  the  same  as  the  air  temperature  in  daytime  mainly  because 
of  its  high  albedo  [%  80%  for  fresh  snow)  and  (3)  responds  rapidly  to 
air  temperature  changes  as  long  as  T\  •  <  0°C,  with  solar  irradiance 

air 

and  wind  speed  exerting  comparatively  smaller  effects.  The  RADC  data 

p 

reflect  (2)  and  (3)  in  that  daytime  values  of  Cn  at  2m  remained  quite 

1 C  0/0 

small.  Typically,  values  were  less  than  10  m'  and  had  small  varia¬ 
bility. 

For  Ta,  >  0°C,  however,  as  long  as  there  is  a  snow  cover  of  sub- 

stantial  depth,  its  surface  temperature  remains  near  0°C  regardless  of 

how  much  warmer  the  air  temperature  becomes.  The  advection  of  T  .jr  >  0°C 

2 

produces  values  of  Cn  that  are  proportional  mainly  to  the  air  tempera¬ 
ture  alone,  with  cloudiness  and  wind  speed  having  secondary  effects.  In 

2 

fact,  one  of  the  largest  increases  in  Cn  in  the  shortest  time  period  for 

the  RADC  data  occurred  at  night  in  an  18-minute  period  between  2254  and 

2 

2312  GMT  on  15  February  1982.  The  3-minute  values  of  Cn  increased  from 
-17  -2/3  -12  -2/3 

near  10  m  to  near  10  m  .  There  were  31  cm  of  snow,  the  2-m 

temperature  was  about  +7°C,  and  the  wind  speed  increased  from  about  0.5 

mph  to  2.1  mph  when  this  occurred. 

It  can  be  expected  that  the  most  common  location  and  time 

for  the  above  condition  to  occur  are  mid-latitude  regions  with  snow  and 

ice  that  is  deep  (thick)  enough  to  last  well  into  springtime.  It  is 

possible,  for  example,  to  have  occasional  periods  with  air  as  warm  as 

about  20°C  moving  over  a  snow  or  ice  surface  until  complete  melting 

occurs.  A  frequent  effect  of  the  cold  surface  on  the  warmer  air  moving 

over  it,  however,  is  to  cause  condensation  as  fog  in  the  first  few 

meters.  Whether  or  not  fog  forms  depends  largely  on  temperature  and 
humidity  characteristics  of  the  air  mass,  but  if  it  does  form,  the 
meaning  and  effects  of  optical  turbulence  are  likely  to  be  obscured. 
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3.2.2  Boundary  layer  profiles 

The  20-minute  averages  of  profile  data  for  8  heights  on  the 

NOAA  BAO  tower  provide  unique  resolution  in  both  time  and  height  re- 

2  2 

garding  the  behavior  of  Cn(Cj )  and  meteorological  variables  between 

10  m  and  300  m.  The  nocturnal  periods  with  a  cloudless  sky  or  scattered 

clouds  listed  in  Table  8  were  used  to  test  the  effectiveness  of  wind 

speed  as  a  variable  to  describe  Cn(Cj)  variations  not  only  at  one  height, 

o 

but  also  vertical  profiles  of  C^,  temperature  and  wind  speed.  Averages 
of  these  variables  for  each  of  the  8  heights  were  calculated  in  succes¬ 
sive  1  ms  ^  categories  of  wind  speed  measured  at  10  m  (u^q)  for  catego¬ 
ries  from  0-1  ms  ^  to  10-11  ms’^ .  The  results  are  shown  i.n  Fig.  13a 

0 

through  13k  for  temperature  and  and  Fig.  14a  through  Fig.  14e  for 

wind  speed.  Temperature  and  wind  speed  plotted  against  height  (z)  are 

2 

shown  with  linear  scales,  and  log  C-p  with  the  same  5-decade  scale  for 
all  figures,  is  shown  against  log  z.  The  category  of  u^Q  is  given  at 
the  top  of  each  figure.  In  the  caption  for  each  figure,  the  range  of 
sample  sizes  for  that  category  of  u-jq  is  given  because  the  sample 
sizes  were  limited  and  varied  not  only  with  height  but  also  from  one 
category  of  wind  speed  to  another.  Table  11  lists  sample  sizes,  maximum 

0  O 

and  minimum  values  of  C^.  at  10  m,  profile  slopes  and  the  top  of  the 
SBL  for  each  u^q  category. 

A  general  feature  of  the  BAO  tower  results  is  that  they  show 
that  changes  in  (1)  actual  values  of  at  10  m,  (2)  slopes  of 
profiles,  (3)  the  top  of  the  SBL  Un),  and  (4)  shapes  and  slopes  of 
temperature  and  wind  profiles  are  affected  by  changes  in  u^q  in  a 
consistent  manner.  For  example: 

1)  c|(10)  increases  from  about  0.03°C  m"^  for  0-1  ms'^  to 


log(CTSQ) 


? 

g.  13a  Temperature  versus  height  (top)  and  log  C,  versus  log 
neight  (bottom)  at  BAO  for  nighttime  condition^  with  a  cloud¬ 
less  sky  or  scattered  cloudiness  and  a  10-m  wind  speed  (u^q) 
category  of  0-1  ms~U  Sample  sizes  range  from  2  to  4. 


log(CTSQ) 


Fit,.  13b 


Same  as  13a  but  with  =  1-? 
sample  sizes  from  13  tou20 * 


-5 


-4 


-3 


-2 


-1 


log(CTSQ) 


Fig.  13c  Same  as  13a  but  with  u,q  =  2-3  ms"^  and  c 
sample  sizes  from  17  tiu 30. 
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2  («)  4  -  5  a/seo 


TEMPERATURE  (C) 


log(Z) 


Fig.  13e  Same  as  13a  but  with  u i q  =  4-5  ms  ^  and  range  in 
sample  sizes  from  34-41. 


lowers*)) 


=  6-7  m$"^  and  a  sample 


Fig.  13g  Same  as  13a  but  with  u,n 
size  of  2.  ,u 


log(CTSQ) 


Fig.  13h  Same  as  13a  but  with  uin 
size  of  2. 


7-8  ms 
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Z  (■)  8-0  «/seo 


lofl(Z) 


TEMPERATURE  (C) 


log(CTSQ) 

Fig.  131  Same  as  13a  but  with  =  8-9  ms’^  and  a  sample  size  of  4. 


log(CTSQ) 


i  Fig.  1 3 j 

k 

♦ 


Same  as  13a  but  with  uln  =  9-10  ms  ^  and 
size  of  4. 


i 


sample 
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2  (■)  10  -  11  b/soo 


log(Z) 


TEMPERATURE  (C) 


log(CTSQ) 


Fig.  13k  Same  as  13a  but  with  uln  -  10-11  ms and  a  sample 
size  of  1.  ,u 


KIND  SPEED  («/sec) 


Fig.  14a  Wind  speed  versus  height  at  BAO  for  nighttime  conditions 
with  a  cloudless  sky  or  scattered  clouds  for  wind  speed 
categories  of  0-1  ms'1  and  1-2  ms-1  corresponding  to 
Fig.  13a  and  b. 


WIND  SPEED  (ft/sec) 


4a  but  for  wind  speed  categories  of  2-3 
is~^  corresponding  to  Fig.  13c  and  d. 


flND  SPEED  (a/sec) 
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Table  11  Cy  maxima  and  minima  at  10m,  profile  slopes, 
and  zn  for  various  uy q  for  the  BAO  profiles 
shown  in  Fig. 13a  through  Fig.  13k. 


u10 
(ms** ) 

Sample  Size 

C2 

(°C 

MAX 

nf2'3) 

MIN 

2 

A  log  Cy 

A  log  z 

zn 

(m) 

0-1 

2-4 

0.03 

0 . 2x1 0~4 

* 

170 

1-2 

13-20 

0.013 

4.3  " 

17u 

2-3 

17-30 

0.012 

5.6  " 

-0.98 

120 

3-4 

13-28 

0.017 

4.1  " 

-1.08 

100 

4-5 

34-41 

0.065 

5.7  " 

-1.67 

60 

5-6 

1 9-2*1 

0.089 

5.5  " 

-1.59 

50 

6-7 

2 

0.156 

3.7  " 

-2.16 

50 

7-8 

2 

0.010 

10  "  (*) 

-0.50  (*) 

50 

8-9 

4 

0.009 

13  "  (*) 

-0.55  (*) 

45 

9-10 

4 

0.007 

15  11  (*) 

-0.45  (*) 

30 

10-11 

1 

0.006 

16  "  (*) 

-0.36  (*) 

25 

not  determinable. 

t*\  ? 

v  '  outlying  small  value  of  Cy  at  200  m 
not  considered. 
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-1  -2/3 

near  0.1  for  u-jq  near  6  ms  and  decreases  to  less  than  0.01  °C  m  ' 

for  lT-j q  >  7  ms~^.  These  changes  are  similar  to  those  for  a  single 
height  discussed  in  Section  3.2.1  in  that  Cy  is  small  at  low  wind  speeds, 
increases  to  a  maximum  as  wind  speed  increases,  and  decreases  above  a 
certain  speed.  The  occurrence  of  a  maximum  near  u^  %  6  ms"^  seems  high 
compared  to  the  .  esults  in  section  3.2.1,  but  it  may  be  due  to  factors 
such  as  height  differences,  spectral  estimates  compared  with  scintillo¬ 
meter  estimates,  and  the  abrupt  decrease  in  sample  size  from  about  20 
samples  for  the  u^  -  5-6  ms'^  category  to  2  for  u^q  ->  6  ms""'. 

2)  Cy  decreases  with  height  for  all  wind  speed  categories, 
with  the  most  marked  decrease  (from  %  10“^°C  m"^  to  ^  10"^°C  m’^) 

occurring  between  100  m  and  300  m  for  u^  =  0-1  ms"^ .  This  height  change 
2 

in  Cy  corresponds  to  the  height  interval  within  which  the  temperature 
profile  changes  from  a  large  inversion  to  nearly  adiabatic.  These 
features,  plus  the  occurrence  of  a  wind  speed  maximum  near  170  in  as 
shown  in  Fig. 14a,  support  the  idea  that  for  u^q  *  0-1  ms”  the  region 
near  170  m  is  the  top  of  the  SBL,  or  zn<  With  increasing  wind  speed, 
the  following  is  indicated  by  the  profiles  regarding  their  slopes,  shapes, 
and  V 

a)  The  negative  slope  of  the  Cy  profiles,  estimated 
for  wind  speeds  beginning  with  u^q  =  2-3  ms"\  increased  negatively 
f**om  about  -1  to  near  -2  for  u y q  =  6-7  ms"^  and  decreased  to  -0.5  or 

less  for  U.JQ  >  7  ms’  .  For  speeds  from  4-6  ms’^  it  was  close  to  -4/3. 

2  -1 
The  decrease  of  Cy  between  100  m  and  300  in  for  u^q  >  6  ms  was  less 

than  half  that  for  u^Q  %  6  ms~*. 
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b)  If  zn  is  defined  by  the  top  of  the  ground-based 
inversion,  a  decrease  from  about  170  m  to  30  m  occurred  as  u-jq  in¬ 
creased  to  %  11  ms  \  contrary  to  eq.  9.  A  similar  lowering  of  the 
height  of  maximum  wind  speed  in  the  3Q0~m  layer  was  indicated. 

c)  A  lack  of  sufficient  samples  for  individual  catego¬ 
ries  for  u-j q  >  5  ms"^  is  a  shortcoming  in  an  interpretation  of  these  re¬ 
sults  but  collectively,  for  u^q  >  6  ms~\the  results  are  consistent. 

The  BAO  results  are  corroborated  by  results  from  the  17 
and  9  CLEAR  I  and  II  soundings,  respectively,  that  were  obtained  in 

O 

stable  and  nearly  cloudless  conditions.  Average  profiles  of  Cp  and 

temperature  from  about  12  m  to  600  m  for  the  CLEAR  I  and  II  soundings 

are  shown  in  Fig.  15  and  Fig.  16,  respectively.  In  each  figure,  the 
2 

coordinates  for  Cn  are  logarithmic  and  those  for  temperature  are  linear. 

For  log  cjj,  the  abscissa  scale  is  the  same  (v  lo"16  to  10~14).  The 
2 

numbers  on  the  Cn  ordinate  are  actual  heights  that  correspond  to  the 
logarithmic  height  scaling  by  computer.  The  average  slope  of  C^  is 
close  to  a  -4/3  slope  that  is  shown  for  reference.  A  dry  adiabatic 
lapse  rate  (r  ^  1°C/I00  m)  is  shown  on  the  temperature  graph  for 
reference.  The  average  height  of  zn,  as  defined  by  changes  in  tempera¬ 
ture  and  relative  humidity,  is  near  120  m  for  CLEAR  I  and  150  m  for 
CLEAR  II. 

2 

for  many  profiles  of  C  ,  a  secondary  maximum  could  be 
identified  near  z  ,  but  its  magnitude  and  vertical  extent  were  smaller 
than  for  a  daytime  CBL.  An  analytical  description  of  it  was  not 
attempted  for  the  SBL.  Instead,  the  -4/3  slope  indicated  by  the 

p  p 

soundings  was  accepted  as  a  reasonable  approximation  of  C^(C^)  through 
the  SBL  up  to  about  2900  m. 


L--V  IT*  L.  \!  ' 
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?  2 

3.3  C  (Ct)  in  near-adiabatic  conditions 
n  T 

For  the  CLEAR  I  data,  the  7  profiles  of  C2  in  Table  4  that 
arc  labeled  "n"  are  those  considered  to  be  representative  of  near¬ 
neutral  ,  or  adiabatic  temperature  profile  conditions  near  the  ground. 

The  profiles  were  measured  near  sunrise  and  sunset  and  are  near  the 

2 

times  when  the  smallest  values  of  C  were  measured  with  the  scintillo- 

n 

meters  at  4  and  14  m. 

2 

An  average  of  the  7  profiles  of  Cp  from  about  12  m  to  2700  m 

is  shown  in  Fig.  17  together  with  profiles  of  temperature  and  relative 

humidity.  ^  -3/4  slope  is  shown  for  reference  on  the  log  Cn  profile  and 

a  dry  adiabatic  lapse  rate  (r)  is  shown  on  the  temperature  profile. 

2  2 

Samnle  sizes  are  also  gm  .  for  Cp.  The  range  of  Cp  values  is  from 

4.2  x  10'16  m  2/3  at  12  m  to  4.2  x  10~18nf2/3  at  about  1500  m.  Fe?- 

2 

tures  of  the  average  Cp  profile  are  (1)  an  apparent  secondary  maximum 

near  2100  m  that  corresponds  to  a  maximum  in  relative  humidity,  with 

some  change  toward  stability  in  thr  temperature  profile  and  (2)  close 

2 

to  a  -4/3  change  of  Cp  with  height  even  though  surface  layer  conditions 
were  nearly  adiabatic.  As  v.'ith  the  temperature  profiles  for  unstable 
conditions,  it  can  be  shown  thai  very  small  changes  ir.  che  profile 
slope  toward  a  more  stable  profile  near  z..  can  produce  significant 
changes  in  the  distribution  and  exchange  of  heat  and  moisture  below 


and  through  z^. 
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4.  Model  testing  and  software  development 

Results  given  in  prior  sections  have  demonstrated  the  effective- 

2 

ness  of  solar  irradiance  and  wind  speed  in  determining  Cn  not  only  at 

2 

a  reference  level  near  the  ground  but  also,  to  some  extent,  Cn  profiles 
through  the  boundary  layer.  Measurements  of  solar  irradiance,  however, 
are  not  routinely  available.  For  the  Cn  model  developed  here,  solar 
irradiance  is  estimated  from  cloud  and  other  information;  cloudiness 
and  wind  speed  are  the  input  variables.  Similarly,  measurements  of 
are  not  routinely  available,  but  estimates  can  be  made  with  cloudiness 
and  wind  speed  information. 

4.1  Unstable  stratification 

The  solar  irradiance  (direct  plus  diffuse  irradiance  on  a 
horizontal  surface)  portion  of  the  model  was  developed  with  the  following 
inputs: 

1)  Latitude  &  longitude, 

2)  Day  of  year, 

3)  Time  of  day, 

4)  Cloud  amount,  preferably  at  the  height  of  each  cloud  layer, 

5)  Cloud  heights 

a)  low-middle  (<  20,000  ft.) 

b)  low-middle  &  high  (mixed  layers) 

c)  high  (>  20,000  ft.). 

Inputs  1)  to  3)  lead  to  calculations  of  hourly  values  of  extra¬ 
terrestrial  solar  irradiance  on  a  horizontal  surface  and  of  the  elevation 


Jk _ ^ -  •* _ 3L 
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angle  of  the  sun.  Inputs  4),  5)  and  solar  elevation  adjust  the  extra¬ 
terrestrial  irradiance  values  for  determining  amounts  reaching  the 
earth's  surface.  Steps  in  the  program  are  as  follows: 

(1)  Extraterrestrial  solar  irradiance  is  calculated  with 
inputs  1),  2)  and  3)  (Paltridge  and  Platt,  1976). 

(2)  An  "equivalent  cloud  amount"  denoted  by  cc  is  determined 
with  inputs  4)  and  5)  as  follows: 

(a)  for  low  to  middle  clouds,  cc  =  total  sky  cover  (in 
tenths) ; 

(b)  for  high  clouds,  cc  =  1/2  total  sky  cover;  and 

(c)  for  mixed  layers:  cc  =  total  sky  cover  -  1/2  high  sky 
cover,  or  cc  =  .75  total  sky  ccver  if  amounts  at 
various  heights  are  not  reported. 

It  is  assumed  that  low  and  middle  clouds  have  similar  diminution  effects 
on  solar  irradiance.  With  an  "equivalent  cloud  amount"  and  the  solar 
elevation  angle,  solar  irradiance  is  determined  with  the  following  empiri¬ 
cal  equations  (Turner  and  Mujahid,  1984),  where  E  is  solar  elevation;  ETRAD 
is  extraterrestrial  irradiance;  and  RAD  is  irradiance  at  the  surface: 

for  E  >  60°;  RAD  «  ETRAD  (.6423  +  .9109(cc)2Sin(E)  -  1.2873(cc)2 

+  1.222  Sin(E)); 

for  E  <20°;  RAD  *  ETRAD  (.308  -  1 .165(cc)2Sin(E)  -  .0586(cc)2 

+  1.0743  Sin(E));  (17) 

for  20°  <  E  <40°:  RAO  *  ETRAD  (.5695  -  .1065(cc)2Sin(E)  -  .4755(cc)2 

+  .2809  Sin(E));  and 

for  40°  <  E  <_60°;  RAD  =  ETRAD  (.7862  -  .2736(cc)2Sin(E)  +  ,6943(cc)2 

-  .0467  Sin(E))  . 
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Values  of  solar  irradiance  (RAD)  calculated  with  the  model 
were  compared  with  measured  values  for  17  days  in  the  CLEAR  I  data.  The 
measured  data  consisted  of  194  values  of  irradiance  on  the  hour  that 
were  abstracted  from  graphs  of  daily  variations  of  irradiance  provided 
for  the  CLEAR  I  experiments.  The  cloud  observations  used  in  the  calcula¬ 
tions  were  available  only  for  Holloman  AFB,  about  68  km  from  where  the 
measurements  were  made.  A  ratio  of  calculated  to  measured  values  (R) 
and  its  inverse  (Inv  R)  were  defined: 


Ratio(R)  = 


RAD  calc 
RAD  meas 


Ir»vR  =  RifR>_l 

»  l  if  R  <  1 


08) 


Table  12  lists  the  individual  days  from  CLEAR  I,  the  number  of  samples 
and  the  average  of  Inv  R  for  each  day  and  Table  13  summarizes  the 
test  results  in  terms  of  R  and  Inv  R.  The  results  given  in  Table  12 
and  13  show  that  the  irradiance  model  underestimates  the  true  irradiance 
more  often  than  it  overestimates  it.  Of  the  194  calculations,  52.6%  of 
them  were  within  90%  of  the  measured  values  and  68.7%  of  them  were 
within  75%  of  measured  values. 

2 

The  same  17  days  were  used  to  test  the  Cn  single-height  model. 
Computations  of  RAD  together  with  wind  speed  are  the  components  of  the 
model  as  given  by  Equation  13  in  Section  3.1.1.  The  measured  data  consisted 
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Table  12  Average  values  of  Inv  R  for  solar  irradiance  for 
CLEAR  I  days 


Date  (1984)  Sample  Size  Inv  R 


8/29 

13 

1.56 

8/30 

8 

2.65 

8/31 

13 

1.14 

9/1 

13 

1.11 

9/7 

13 

1.13 

9/8 

13 

1.16 

9/10 

13 

1 .34 

9/11 

12 

1.66 

9/18 

12 

1.43 

9/20 

12 

1.10 

9/21 

12 

1.23 

9/22 

8 

1.57 

9/23 

12 

1.22 

9/24 

12 

1.30 

9/25 

11 

1.97 

9/26 

6 

1.52 

9/27 

11 

6.78 

Tablel3  Results  of  testing  the  solar  irradiance  model 

against  measurements  for  CLEAR  I  for  194  samples. 
R  is  the  ratio  of  calculated  to  measured  values 
of  solar  irradiance 


Min. 

Max. 

Avq- 

Sample  Size 

%  of  Samples 

R  <  1 

.02 

.99 

.85 

106 

54.6% 

R  *  1 

1.00 

1.00 

1.00 

5 

2.6% 

R  >  1 

1.01 

8.84 

1.55 

83 

42.8?, 
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of  tower  1,  HIDL  values  on  the  hour  of  (1)  wind  speed  at  4  m  (2)  solar 

2 

irradiance  and  (3)  C  for  the  4-m  scintillometer.  For  the  calculations 

n 

2 

of  Cn>  RAD  was  estimated  as  described  above  and  a  10-m  wind  speed  ( WS-j q) 
was  estimated  with  the  4-m  speed  (W$4)  with  the  power  law  relationship: 

WS1Q  S  WS4  {^)0'209  S  1.211  WS4  ,  (19) 

where  the  power  law  exponent  for  unstable  conditions  is  an  interpolated 

value  based  on  work  by  Touma  (1977).  This  approximation  for  WS-jq  is 

considered  sufficient  in  view  of  the  comparatively  much  larger  effect 

of  RAD  on  C2. 

n 

2 

The  ratios  of  calculated  to  measured  values  of  C  are  as  de- 

n 

fined  in  Eq.  18.  Table  14  lists  the  individual  days,  sample  sizes  and  the 

2 

average  of  Inv  R  for  Cn  for  each  day  and  Table  15  summarizes  the  results. 

2 

They  show  that  the  model  overestimates  nearly  as  often  as  it  under¬ 
estimates  it.  In  addition,  they  show  that  67.9%  of  the  estimated  values 

are  within  a  factor  of  2  of  the  measured  values.  As  would  be  expected, 

2 

the  best  agreement  for  both  RAD  and  Cn  was  found  for  periods  with  no 
cloudiness. 

2 

The  entire  Cn  model  for  unstable  stratification,  from  estimates 

2 

of  Cn  for  a  single  height  with  cloudiness  and  wind  speed  information  to 

estimates  of  boundary  layer  profiles  with  the  K-T  model  and  estimates  of 

2 

z..  was  tested  against  3  measured  Cn  profiles  for  the  CLEAR  I  period  and 
3  for  the  CLEAR  II  period.  Table  16  lists  the  dates  ano  times  of  the 
soundings  along  with  relevant  information  defined  below.  The  symbols  are 
as  used  in  the  software  for  the  model. 
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Table  14  Average  values  of  Inv  R  for  for  CLEAR  I  days 


Date  (1984) 

Sample  Size 

Inv  R 

8/29 

13 

3.66 

8/30 

7 

3.53 

8/31 

13 

1.86 

9/1 

13 

1.82 

9/7 

11 

1.22 

9/8 

11 

1.49 

9/10 

12 

2.68 

9/11 

10 

2.19 

9/18 

13 

2.52 

9/20 

13 

2.20 

9/21 

12 

1.68 

9/22 

7 

3.09 

9/23 

13 

3.22 

9/24 

13 

2.56 

9/25 

9 

10.47 

9/26 

13 

21.55 

9/27 

13 

8.32 

2 

Table  15  Results  of  testing  the  single-height  Cn 
model.  R  is  the  ratio  of  calculated 
with  irradiances  and  wind  speeds  to 
measured 


Min. 

Max. 

Avg. 

Sample  Size 

%  of  Samples 

R  <  1 

.10 

.98 

.64 

105 

53 . 6ii 

R  =  1 

1.00 

1.00 

1 .00 

1 

.5*5 

R  >  1 

1.01 

98.06 

7.18 

90 

45. 9--; 

R  <  .5 

.10 

.45 

.25 

17 

8.7S 

R  <  2 

.50 

1.98 

.94 

133 

67.9% 

R  >  2 

2.02 

98.06 

12.72 

46 

23 . 5S 
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Table  16  Information  relevant  to  CLEAR  I  and  II 
soundings  to  test  the  profile  model. 


DATE  (1984) 

TIME 

CC 

CH 

WS 

(est) 

(mea) 

840907 

1417 

0 

0 

3.9 

1180 

2100 

840908 

1343 

1 

5 

1.8 

950 

2600 

840910 

1457 

3 

1 

3.5 

1080 

2700 

850227 

0938 

2 

7 

3.6 

820 

450/1000 

850301 

0950 

1 

1 

2.0 

580 

1700 

850305 

1441 

<1 

7 

5.1 

1450 

1100 

Time  is  mountain  standard  time. 

CC  is  cloud  amount  in  tenths  of  opaque  cloudiness, 

CH  is  the  cloud  height  code  given  in  the  radiance 
model, 

WS  is  wind  speed  (ms-1)  at  10m  for  CLEAR  I  and  at 
4m  for  CLEAR  II,  and 

z-  is  the  mixed  layer  height  (m)  estimated  from 
'  the  Smith  (1977)  nomogram  (est)  and  from  tempera¬ 
ture  and  humidity  profiles  (mea). 

2 

The  model  provided  estimates  of  Cn  at  4m,  with  values  at  10,  20, 

40,  60m  ...  etc.  calculated  with  eq.  14.  The  calculations  were  then 
averaged  in  the  same  manner  as  the  measured  soundings.  Results  are  shown 
in  Fig.  18a  through  18f.  The  solid  line  is  the  model  prediction  and  the 
dots  are  measured  values.  Except  for  the  profile  on  9/10/84,  where 
cloudiness  estimates  may  have  been  erroneous,  reasonable  agreement  between 
measured  and  estimated  values  can  be  noted.  A  major  source  of  error  is 
the  large  discrepancies  between  estimated  and  measured  values  of  z.  shown 
in  Table  16.  Reasons  for  the  differences  were  not  established  as  part  of 
the  present  work. 
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? 

Fig.  18a  Comparison  of  measured  profiles  of  C  in  unstable 

stratification  (CLEAR  I  and  II)  witnprofiles  calcu¬ 
lated  from  cloudiness,  wind  speed,  location,  time 
and  date  information  for  1417  MST,  9/7/84  and 
1343  MST,  9/8/84. 


Fig.  18b  Same  as  18a,  but  for  1457  MST,  9/10/84  and 
1438  MST,  2/27/85. 
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4.2  Stable  stratification 

2 

As  discussed  in  Section  3.2.2,  the  decrease  in  Cy  with  height 

appeared  to  be  inversely  dependent  on  u^g,  but  there  were  questions  for 

—  -1  2  —  -1 
u-jg  <  3  ms  because  of  curvature  in  the  Cy  profile  and  for  i.:-jg  >  6  ms 

because  of  small  sample  sizes.  For  4  <  u1ri  <  6  ms-1,  where  a  fully 

turbulent  regime  is  likely,  however,  the  slope  was  near  -4/3.  The  average 

of  the  soundings  for  CLEAR  I  and  II  also  indicated  a  slope  near  -4/3  above 

the  -Mrst  sounding  height. 

Based  on  these  experimental  results,  a  -4/3  slope  was  accepted 

2  2 

in  the  present  work  for  describing  Cn(Cy)  in  a  turbulent  stable  boundary 
layer  above  10m.  A  -2/3  dependence  predicted  by  similarity  theory  may 
exist,  but  as  indicated  by  Walters  and  Kunkel  (1981),  it  is  probably 
below  the  lowest  measurement  height  which,  in  their  case,  was  8  m. 

The  final  model  for  stable  stratification,  then,  consists  of 
the  Cn~wind  speed  regression  equation  obtained  for  the  BAO  data  for  10  m 
coupled  with  a  C2n  *  z~^3  variation  above  10  m.  The  model  was  tested 
with  3  CLEAR  I  and  3  CLEAR  II  soundings  selected  arbitrarily  and  the 
results  are  shown  in  Fig.  19a  through  19c  for  heights  from  10  m  to  3000  m. 
The  solid  line  is  the  model  prediction. 

It  can  be  noted  that  for  most  cases,  below  about  800  m, estimated 
values  are  greater  than  measured  values,  but  a  -4/3  slope  provides  a 
reasonable  estimate  of  C^(z).  There  could  be  many  reasons  for  the  be¬ 
havior  and  differences  shown.  They  involve  terrain  effects  that  are 
especially  important  in  stable  conditions,  cloudiness  differences  between 
Holloman  AFB  and  the  measurement  site  (for  the  CLEAR  1  data)  and  other 
effects  that  are  less  evident  but  possibly  as  important. 


-17.4  -1G3  - 15.2  -14.1  -110 
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Fig.  19a  Comparison  of  measured  profiles  of  Cp  in  stable 

stratification  (CLEAR  I  and  II)  with  profiles  calcu¬ 
lated  from  the  BAO  C^-wind  speed  relationship  and  a 
-4/3  height  variation;  2113  MST,  9/3/84  and  2056  MST 
9/6/84. 
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Fig.  19b  Same  as  19a,  but  for  2005  MST,  9/10/84  and 
1815  MST,  2/28/85. 


asms* 
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4.3  Fortran  Program 

CNSQ  PROG  and  STABLE  are  FORTRAN  IV  programs  that  were 

p 

developed  to  (1)  calculate  Cn  at  one  height,  and  with  additional 
information,  (2)  calculate  its  height  variation  for  the  boundary 
layer.  The  programs  are  described,  together  with  a  discussion  of  input 
variables,  operation  and  output  variables,  and  adaptation  for  use  at 
different  sites.  Examples  of  how  they  are  run  on  the  FORTRAN  IV 
compiler  on  the  Michigan  Terminal  System  are  also  given  with  soundings 
for  White  Sands,  NM*  listed  in  Appendix  A  and  B. 

4.3.1  Input  variables  for  unstable  stratification 

The  following  are  the  input  variables  for  the  program 
that  are  necessary  for  its  operation: 

a)  Year,  month,  day,  and  time, 

b)  Cloud  cover  and  cloud  height  code, 

c)  Wind  speed  at  10  meters  and 

d)  Height  of  the  mixed  layer  (z^). 

a)  Year,  month,  day,  and  time.  Because  the  program  is  pre¬ 
sently  written  for  White  Sands,  NM,  time  is  in  Mountain  Standard  Time. 
The  program  listed  is  valid  for  one  hour  after  sunrise  to  one  hour  be¬ 
fore  sunset. 

b)  Cloud  cover  and  cloud  height  code.  Cloud  cover  is  in 
tenths  of  total  sky  cover  (preferably  opaque  sky  cover)  and  the  cloud 
height  codes  are  described  below: 

0  -  clear 

1  -  low  clouds;  cloud  base  <  10,000  ft 

2  -  middle  clouds;  10,000  ft  <_  cloud  base  <  20,000  ft 

3  -  low  &  middle  clouds 


4  -  low  &  high  clouds 

5  -  low,  middle  &  high  clouds 

6  -  middle  &  high  clouds 

7  -  high  clouds;  cloud  base  20,000  ft 

c)  Wind  speed  at  10  meters.  Because  wind  speed  is  routinely 

reported  in  knots  as  a  measured  variable  in  conventional  surface  weather 

2 

observations,  this  unit  was  used  in  the  program  for  C  .  If  wind  speed 
is  available  for  a  height  other  than  10  m,  the  following  form  of  the 
power  law  for  a  moderately  unstable  boundary  layer  can  be  used  (Touma,  1977) 

WS1Q  £  WS£ ( 1 0/z ) 0  * 209  , 

where  z  is  the  height  of  the  measured  wind  speed  in  meters. 

d)  Height  of  the  mixed  layer  (z^.  The  time  of  day,  cloud 
cover,  and  average  wind  speed  at  10  m  (or  an  estimate  of  it,  as  above) 
are  required  for  z-  estimates  with  the  Smith  (1977)  nomogram.  The 
mixed  layer  height  estimated  for  the  example  given  here  to  illustrate 
the  operation  of  the  program  was  about  950  m. 

4. 3. 1.1  Description  and  operation 

The  following  is  a  brief  description  of  the  variables, 
statements, and  functions  in  CNSQ.PR0G: 

Lines  1-8;  RLAT :  latitude  of  site  in  0 
RL0N:  longitude  of  site  in  0 

RZONE:  15°  longitude  increment  for  the  time  zone  of  the  site 
(105°  for  MST) 

SCON:  The  solar  constant  in  ly/min. 

EXPO:  -4/3  (the  slope  of  the  variation  of  with  height) 

2 

H4:  4  meters,  the  height  for  which  a  "surface  Cn"  value 
will  be  estimated. 
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Lines  9-24:  enter  input  variables. 

Line  25:  IDATE  is  the  date  in  16  format  ( YR/MON/DAY) 

Line  26-41:  Calculate  the  zenith  angle  and  elevation  angle  of  the 
sun  for  the  time  desired,  with  which  a  value  of  extra¬ 
terrestrial  solar  irradiance  is  calculated. 

RTIME:  time  in  hours  &  tenths  of  an  hour.  For  example 
1344  =  13.7333. 

EQTM:  the  value  from  the  "equation  of  time",  which  accounts 
for  the  change  in  solar  noon  as  a  function  of  time 
of  year. 

STIME:  solar  time. 

HAS:  hour  angle  of  the  sun  in  radians. 

DECL:  solar  declination  in  radians. 

COSZ:  cosine  of  the  solar  zenith  angle. 

ETRAD:  solar  irradiance  at  the  "top  of  the  atmosphere"  In 
ly/min. 

Z:  solar  zenith  angle  in  degrees. 

E:  solar  elevation  angle  in  degrees. 

Lines  42-60:  Calculate  global  solar  irradiance  at  the  earth's  surface 
with  cloud  cover,  cloud  height  code,  and  solar  elevation 
angle  as  inputs. 

RADC:  calculated  global  solar  irradiance  at  the  earth's 
surface  in  ly/min. 

Line  61:  Converts  wind  speed  from  knots  to  m/sec. 

Line  62-63:  The  regression  eqn.  for  estimating  Cn  from  calculated  ra¬ 
diance  and  wind  speed. 

-8  -1/3 

Line  64:  sets  Cn  =  10  m  '  if  the  regression  equation  produces  a  value 
below  10’8. 


v»  a  v*  WkTik'X  V?  VV  -V  A  vV.V:\W.  .-A  •■V--v 
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Line  65:  Calculates  Cn  where  CNSQC  is  calculated  Cn  in  m  7  at  4  m. 

2 

Lines  66-67:  Write  the  "surface  measurement"  of  C  ,  along  with  the 
input  variables,  into  a  file  called  "SURFACE"  in  the 
sample  run. 

2 

Lines  68-77:  Compute  the  Cn  profile  with  equation  14.  Lines  72-73 

2 

calculate  1 og^ Q  of  Cn  and  height,  and  line  76  adjusts 

.the  height  interval  (20  m  in  this  case).  Line  69  deter- 

2 

mines  how  many  values  of  Cn  will  be  calculated.  In  the 

example,  there  are  151  values  from  20  m  to  3000  m. 

HGT:  height  in  meters. 

CNZC:  calculated  C2  in  m"2/3. 

n 

HGTL0G:  1og1Q  HGT. 

CNZL0G:  log1Q  C2. 

Lines  81-110:  A  subroutine  used  in  line  26  which  converts  the  year,  month, 

and  day  into  Julian  calendar  day,  which  is  then  used  to 

compute  solar  declination. 

IB:  Julian  calendar  day. 

4. 3. 1.2  Output  variables  and  program  adaptations 

The  output  of  the  program  is  as  shown  for  the  sample  run, 

2  2 

with  values  of  height,  Cn,  log^g  (height),  and  log^  Cp  printed  out.  The 

log^Q  values  are  retained  for  model  output  graphs  of  values  of  log^Q  (height) 
2 

vs.  log^g  C  •  Some  of  the  statements  in  the  program  must  be  altered  to  be 
applicable  to  a  site  other  than  White  Sands,  NM.  A  list  of  these  state¬ 
ments  is  given  below. 

Line  1 :  RLAT 
Line  2:  RL0N 

Line  3:  RZ0NE,  the  15°  longitude  line  corresponding  to  the  time  zone  of 
the  site.  For  example,  for  the  Eastern  Time  Zone,  RZ0NE  s  75. 


k! 
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In  addition; 

Line  69  &  76;  The  do-loop  parameter  and  height  interval  can  be  adjusted 
as  desired. 

4.3.2  Input  variables  for  stable  stratification 

The  FORTRAN  program  CNSQ.  STABLE  was  developed  to  cal- 

2 

culate  the  profile  of  Cn  for  stable  conditions.  It  is  a  relatively  short 

program  and  uses  the  "normal  curve"  approximation  for  the  variation  of 
o 

Cn  with  wind  speed  at  10  meters  with  eq.  (16)  developed  for  the  BAO  data. 

If  wind  speed  is  measured  at  a  height  other  than  10  meters,  the  following 
equation  can  be  used  to  obtain  the  10-m  speed  for  stable  conditions 
(Touma,  1977): 

wslo  =  wsz  (y)'414  ’  (20) 

where  WS-jq  is  the  estimated  wind  speed  at  10  m  and  WSz  is  the  speed 
measured  at  height  z. 

4.3.2. 1  Description  and  operation 

A  brief  description  of  the  statements  in  CNSQ.  STABLE 

is  given  below: 

Line  1:  the  exponent  (-4/3)  used  in  the  equation  for  the  variation  of 
with  height. 

Lines  2-5:  statements  that  enable  the  user  to  enter  the  wind  speed  at  10  m. 
Line  6:  a  conversion  for  wind  speed  from  knots  to  ms~^. 

o 

Line  7:  calculates  an  estimate  of  Cn  at  4  m  with  equation  (16)  where 
CN  SQC  is  C^(4m) . 

2 

Lines  8-9:  write  the  estimated  "surface  C  ". 

n 

2 

Lines  10-18:  involve  a  do-loop  that  calculates  the  Cn  profile.  The  do- 
loop  parameter  (I  in  Line  11)  and  the  height  interval  (20  m 


in  Line  17)  can  be  changed  as  desired.  Lines  15-16  write 
the  following  variables  into  a  file,  where 
HGT  =  height 

CN2C  =  calculated  cf  at  HGT 
HGTLOG  =  log10(HGT) 

CN2L0G  =  log10(C|;). 

2 

In  the  example,  the  "surface"  estimate  of  Cn  is  written  into  the  file 
SURFACE,  and  the  profile  is  written  into  the  file  PROFILE.  The  example 
used  corresponds  to  the  White  Sands,  NM  thermosounding  at  21:13:31  on 
9/3/84. 
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5.  Conclusions  and  recommendations 

1)  For  daytime  conditions  with  half  or  less  of  the  sky  cloud-covered, 

? 

Cn  at  one  height  was  describable  in  relation  to  solar  irradiance  and  wind 
2 

speed.  Cn  increased  with  solar  irradiance  and  wind  speed  until  certain 

values  were  reached  and  was  expressible  in  terms  of  a  third-order  poly- 
2 

nomial  involving  C  and  these  two  variables. 

3  n 

2 

2)  The  Kukharets  and  Tsvang  profile  model  for  Cp  tested  satisfactorily 

2 

against  thermosonde  measurements  of  Cn  through  a  convective  boundary  layer 

up  to  3000  m  provided  that  values  of  mixing  depth  were  obtainable  to  the 

2 

nearest  100  m  and  that  Cn  at  a  reference  height  was  known. 

3)  By  applying  a  model  for  calculating  solar  irradiance  with  cloud 

2 

and  other  information,  and  by  using  the  value  of  Cn  calculated  with  the 

solar  irradiance-wind  speed  regression  equation,  it  was  possible  to 

2 

estimate  Cp  at  a  reference  height.  With  this  information,  together  with 

estimates  of  mixing  height  and  the  Kukharets  and  Tsvang  model,  a  goal 

2 

of  describing  the  variation  of  Cn  through  the  atmospheric  boundary  layer 
with  conventional  meteorological  variables  was  reached  for  unstable 
stratification.  A  major  source  of  error  was  in  estimating  z^. 

4)  For  nighttime  conditions  with  half  or  less  of  the  sky  cloud- 

2 

covered,  Cn  at  a  reference  height  was  describable  in  terms  of  wind  speed 
alone.  Actual  values  depend  on  the  reference  height,  but  in  general, 

Cn  increased  sharply  as  wind  speed  increased  from  nearly  calm,  it 
reached  maximum  values  in  a  narrow  wind  speed  interval,  and  decreased 
at  higher  speeds.  These  changes  were  expressible  with  appropriate 

coefficients  in  the  equation  for  a  normal  curve. 

2  2 

5)  The  Cn-wind  speed  equation  for  calculating  Cn  at  a  reference 


4  w  i  ^  ^  ~  >■>  w  ,»> 
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o 

height,  together  with  a  -4/3  height  variation  of  Cn  above  the  reference 

height  obtained  from  thermosonde  data,  enabled  a  goal  of  describing  the 
2 

variation  of  Cn  through  the  atmospheric  boundary  layer  with  conventional 
meteorological  variables  to  be  reached  for  stable  stratification. 

O 

6)  For  near  adiabatic  conditions,  for  a  reference  height  in  the 

surface  layer  was  at  minimum  values.  Above  about  10  meters,  a  -4/3  height 

variation  through  the  boundary  layer  was  indicated  even  though  very  small 
2 

values  of  C  were  observed, 
n 

p 

7)  Occasionally,  large  morning-to-afternoon  differences  in  for 
the  same  value  of  solar  irradiance  were  observed  mainly  for  the  RADC  data 
but  also,  to  a  lesser  extent, for  other  locations.  They  were  not  fully 
explainable  because  the  necessary  types  of  measurements  were  not  availa¬ 
ble.  This  characteristic  is  probably  the  major  source  of  error  in  the 

2 

Cn-solar  irradiance-wind  speed  relationship.  For  conditions  over  sur¬ 
faces  with  high  soil  moisture  and  large  evaporation,  two  relationships, 
one  for  morning  and  one  for  afternoon,  would  probably  be  more  represen¬ 
tative.  It  is  recommended  that  this  important  question  be  answered  with 
measurements  of  optical  turbulence  over  moist  (grass-covered)  surfaces 
that  are  documented  with  measurements  of  temperature,  humidity  and 
wind  speed  profiles  as  well  as  soil  moisture  and  radiation. 

8)  Significantly  different  values  of  C^(Cy)  are  obtainable  with 
different  methods  of  measurement .  Not  discussed  here  are  yet  other 
methods,  such  as  SODAR  and  LIDAR  that  are  being  used  for  this  purpose. 

It  is  recommended  that  an  intercomparison  experiment  be  conducted  with 
the  various  methods  to  determine,  for  various  turbulence  conditions, 
nOv  only  how  they  compare  but  also  which  method  provides  the  most  re¬ 
presentative  values  of  optical  turbulence  as  it  affecis  the  performance 
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of  optical  systems. 

9)  The  thermosonde  data  used  as  a  reference  in  comparing  profile 
2 

measurements  of  Cn  with  model  estimates  were  invaluable,  but  one  short¬ 
coming  that  limited  their  usefulness  was  a  lack  of  a  standard  lowest 
height  of  measurement.  It  is  recommended  that,  if  possible,  a  thermo¬ 
sonde  system  equivalent  to  that  being  launched,  together  with  a  wind 
speed  and  direction  system,  be  operated  at  one  height  within  the  first 

few  meters  before  and  during  each  thermosounding  to  provide  a  reference 
2 

measurement  of  Cn  and  wind  velocity. 
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Appendix  A.  Fortran  IV  program  for  Cp  model  for  unstable 
stratification,  with  example  for  1344  MST* 

8  Sept.  1984  for  White  Sands,  NM. 


CNSQ.PROG 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 


RLAT-32.4 
RL0N-1O6.37 
RZ0NE-105, . 

PI*3. 141592654 
DTR-PI/180. 

SCON= 1 370 » /697 «  8 
EXP0=-4./3. 

H4«4. 

URITE<  6,101)*' 

101  FORMAT ( * tEnter  wear,  aonth,  daw,  and  tiae  (MST>  <312, 15):  ') 

READ<5, 102)  IYR»IMON,IDAY*ITIME 

102  FORMAT (3X2, 15) 

URITE(6,201> 

201  FORMAT < 'tEnter  cloud  cover  (in  tenths)  and  cloud  height  code  (212)!*) 
READ (5, 202)  XCC,ICH 

202  F0RMAT(2I2) 

WRITE(6,301 ) 

301  FORMAT  ( *  tEnter  wind  speed  at  10  aetres  (in  knots)  (F4.D1  * ) 

READ (5* 302)  US 

302  F0RMAT(F4.1) 

URXTE(6,401) 

401  FORMAT ( ' tEnter  estiaated  height  of  nixed  lauer  (in  aetres)  (F5.01J #) 
READ(5,402)  ZI 

402  FORMAT (F5.0) 

IDATE-I YR* 1 0000+ IMON* 1 00+ ID AY 
CALL  IYRC(ID»IYR,IMONf IDAY) 

RMIN-ITIME-ITIME/1004100 

RTIME-ITIME/100+RMIN/60. 

EQTH»-2.721175+IF*(-.5070817+ID*(.50B2854E-02+I8*(-.2113816E-04 

1  +IB*< .1543856E-05+IB*(-.2542001E-07+IB*( ,1692641E-09+IB*(- 

2  «5555747E-12+IBt ( « 8987469E-15+IBJM-.  57443 1 8E-10)  )  )  )))))) 
STIHE*RTIME-( (RLON-RZONE )/15« -E0TM/60. ) 

HAS-(STIME-12.)*15.*DTR 

DSs( IB-1 )t2. tPI/365. 

DECL«.OO6910-,399912*COS<OS>+.O7O257*SIN<DS)-.OO6758*COS<2.*DS) 

1  +.000907*SIN<2.*DS)-.002697*COS(3.*DS)+.00148*SIN<3,*DS) 

RLAT-RLATSBTR 

COSZsSIN<RLAT)SSIN(DECL)+COS(RLAT)*COS(DECL  )*COS(HAS) 

ETRAD*SC0N*C0SZ 
Z»ARCOS ( COSZ > /DTR 
E“90.-Z 

IF(ICH  *GE.  4  .AND.  XCH  .LE.  6)  GO  TO  24 
IF ( I CM  ,E0.  7)  GO  TO  25 
CC«1.«ICC 
GO  TO  26 

24  CC-.754ICC 
GO  TO  26 

25  CC= . 541CC 
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49 

50 

51 

52 
|  53 

54 

55 

56 

57 
5JB 

59 

60 


I' 6 


1 


6\ 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

71 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 


fl  05 


27 


28 


29 

31 


CC=.1*CC 
IF<E  .LE.  20.)  GO  TO  27 
IF<E  .LE.  40.)  GO  TO  28 
IF(E  .LE.  60.)  GO  TO  29 

RADC«ETRAD*<  ,64234. 9109*CC**2*SIN<E*DTR>-1.2873*C£**2+.1222*SIN<E*DTR) 
GO  TO  31 

RADC-ETRAD*  < . 308-1 ♦ 165*CC**2#SIN<  E*DTR ) - . 0586*CC**2+1 . 0743*SIN  <  E*DTR ) 
GO  TO  31 

RADC«ETRAD*<.5695-.1065*CC**2*SIN<E*DTR>-»4755*CC**2+.2809*SIN<E*DTR> 
GO  TO  31 

RADC-ETRAD*  <  ,  7862-  ♦  2736*CC**2*SI  N  <  EtDTR  >  +-.6943*CC**2-  .  0467*SIN  < E*DTR ) 
IF<RADC  .LT.  0.)  RADC-O. 


100 


200 


10 


US=US*. 51479 

CNC=«  12724E-06+.  17314E-06*RADC-  .33899E-07*l«-'S+«6223bE-06*RADC**2 
+.10264E-07*WS**2-.39727E-06*RADC**3-.53824E-09*US**3 
IF<CNC  .LT.  l.E-08)  CNC-l.E-08 
CNSQC=CNC**2 

UR I TE <  7 » 1 00  f  I DATE » I T IME » I CC  r I CH » US » RADC » CNSQC 
FORMAT  <I6»I5»2I4»F7.1fF9.2»E14«5) 

HGT-10. 

DO  10  1=1 t 151 

CN2C-CNSQC*  < .046* (HGT/2I )**EXP0+.6*EXP( -' 2*(HGT/ZI-1 . 1 )**2) ) 

/ <  «046*(H4/2I )**EXPO) 

HGTL0G=AL0G10(HGT ) 

CN2L0G=AL0G10<CN2C ) 

WRI TE <  8 . 200 )  HGT »  CN2C » HGTLOG . CN2L0G 
FORMAT  (F5.0»E15.6»3X»2F12.5> 

HGTstI*20. 

CONTINUE 

STOP 

END 


r> 


8 


10 


11 


12 

13 

14 


SUBROUTINE  I YRC< IB» IYR* I MON. I DAY ) 
GO  T0(l»2»3f4»5»6»7»8»9»10»llrt2) 
IB=0 

GO  TO  14 
IB=31 
GO  TO  14 
I  B=59 
GO  TO  13 
I  8=90 
GO  TO  13 
I D=120 
GO  TO  13 
18=131 
GO  TO  13 
18=181 
GO  TO  13 
18=212 
GO  TO  13 
I 8=243 
GO  TO  13 
18=273 
GO  TO  13 
I B=304 
GO  TO  13 
1 8  =  334 

IF<M0D(1YR*4).E0.0>  18=18+1 
18=18+ IDA Y 

IF(IDAY.EO.O)  18=18+15 

RETURN 

END 


I  MON 


EXAMPLE 


♦RUN  *FTN  SCARDS=CNSQ.PROG 
♦Execution  besins  14136105 
No  errors  in  MAIN 
No  errors  in  2YRC 

♦Execution  terminated  14236*07  T*0.148  40*09 

♦RUN  -LOAD  7*SURFACE  8«PR0FILE 
♦Execution  begins  14236224 

Enter  veer*  month*  day*  and  time  (MST)  (3I2»I5)1  840908  1344 

Enter  cloud  cover  (in  tenths)  and  cloud  heidht  code  (212)2  1  8 

Enter  wind  s^-eed  at  10  metres  (in  knots)  (F4«l)2  1*8 

Enter  estimated  heiaht  of  mixed  layer  (in  metres)  (FS*0)2  947* 

♦Execution  terminated  14236244  T*0«054  40*04 

♦ 


♦LIST  SURFACE 

1  840908  1344  1  5  0*9  1*20  0.27228E-12 

♦ 


♦LIST  PROFILE 

z(m) 

c*(»-y3) 

log10(z) 

,o9)0  i 

1 

10. 

0.802465E-13 

1.00000 

“13.09557 

2 

20. 

0.318458E-13 

1.30103 

-13.49695 

3 

40. 

0. 1263B0E-13 

1.60206 

-13.89832 

4 

60. 

0.736023E-14 

1.77815 

-14.13311 

S 

80. 

0. 501541E-14 

.  1.90309 

-14.29969 

6 

100. 

0.372473E-14 

2.00000 

-14.42890 

7 

120. 

0.292094E-14 

2.07918 

-14.53448 

8 

140. 

0.237828E-14 

2.14613 

-14.62374 

9 

160. 

0. 199044E-14 

2.20412 

-14.70105 

10 

180* 

0. 1701 22E-14 

2.25527 

-14.76924 

11 

200. 

0. 147834E-14 

2.30103 

-14.83023 

12 

220. 

0.130204E-14 

2.34242 

-14.88538 

13 

240. 

0.115961E-14 

2.38021 

-14.93569 

14 

260. 

0.104251E-I4 

2.41497 

-14.98192 

15 

280. 

0.944834E-15 

2.44716 

-15.02464 

16 

300. 

0.862398E-15 

2.47712 

-15.06429 

17 

320* 

0.792157E-15 

2.50515 

-15.10119 

18 

340. 

O.731B70E-15 

2.53148 

-15.13556 

19 

360. 

0.679914E-15 

2.55630 

-15.16755 

20 

380* 

0.635051E-15 

2.57978 

-15.19719 

21 

400. 

0.596418E-15 

2.60206 

-15.22445 

22 

420. 

0.563419E-15 

2.62325 

-15.24917 

23 

440. 

0.535691E-15 

2.64345 

-15.27108 

24 

460. 

0.513075E-15 

2.66276 

-15.28982 

25 

480. 

0.495595E-15 

2.68124 

-15.30487 

26 

500. 

0, 483449E-15 

2.69897 

-15.31565 

27 

520. 

0.476989E-15 

2.71600 

-15.32149 

-  131 


26 

540. 

0 « 476718E-15 

2.73239 

>15.32174 

29 

560  * 

0.4B3264E-15 

2.74819 

-15.31581 

30 

580. 

0.497366E-15 

2.76343 

-15.30332 

31 

600. 

0.519841E-15 

2.77815 

-15.28413 

32  • 

620. 

0.551550E-15 

2.79239 

-15.25841 

33 

640. 

0.593347E-15 

2.80618 

-15.22669 

34 

660. 

0.646028E-15 

2.81954 

-15.18974 

35 

680. 

0.710264E-15 

2.83251 

-15.14858 

36 

700. 

0.786531E-15  - 

2.84510 

-15.10429 

37 

720. 

0.875037E-15 

2 .85733 

-15.05797 

38 

740. 

0.975655E-15 

2.86923 

-15.01070 

39 

760. 

0. 108785E-14 

2.88081 

-14.96343 

40 

780. 

0.121062E-14 

2.89209 

-14.91699 

41 

800. 

0.134249E-14 

2.90309 

-14.87209 

42 

.820  . 

0. 148147E-14 

2.91381 

-14.82931 

43 

840. 

0. 162506E-14 

2.92428 

-14.78913 

44 

860. 

0. 177034E-14 

2.93450 

-14.75194 

45 

880. 

0.191402E-14 

2.94448 

-14.71805 

46 

900. 

0.205251E-14 

2.95424 

-14.68771 

47 

920. 

0. 218214E-14 

2.96379 

-14.66112 

48 

940. 

0.229921E-14 

2.97313 

-14.63842 

49 

960. 

0.240023E-14 

2.98227 

-14.61975 

50 

980. 

0.248204E-14 

2.99123 

-14.60519 

51 

1000. 

0. 254198E-14 

3.00000 

-14.59483 

52 

1020. 

0.257801E-14 

3*00860 

-14.58871 

53 

1040. 

0 . 258B84E-14 

3.01703 

-14.58689 

54 

1060. 

0. 257398E-14 

3.02531 

-14.58940 

55 

1080. 

0.253375E-14 

3.03342 

-14.59623 

56 

1100. 

0.246930E-14 

3.04139 

-14.60743 

57 

1120. 

0. 238251E-14 

3.04922 

-14.62296 

58 

1140. 

0. 227593E-14 

3.05690 

-14.64284 

59 

1160. 

0. 215260E-14 

3.06446 

-14.66704 

60 

1180. 

0.201596E-14 

3.07188 

-14.69552 

61 

1200. 

0.186964E-14 

3.07918 

-14.72824 

62 

1220. 

0.171731E-14 

3*08636 

-14.76515 

63 

1240. 

0.156254E-14 

3.09342 

-14.80617 

64 

1260. 

0. 140864E-14 

3.10037 

-14.85119 

65 

1280. 

0. 125859E-14 

3.10721 

-14.90011 

66 

1300. 

0. 111490E-14 

3.11394 

-14.95276 

67 

1320. 

0 . 979605E-1 5 

3.12057 

-15.00895 

68 

1340. 

0.854234E-15 

3.12710 

-15.06842 

69 

1360. 

0.739805E-15 

3.13354 

-15.13088 

70 

1380. 

0.636869E-X5 

3.1398B 

-15.19595 

71 

1400. 

0 . 545549E-15 

3.14613 

-15.26316 

72 

1420. 

0. 465626E-15 

3.15229 

-15.33196 

73 

1440. 

0. 396575E-15 

3.15836 

-15.40168 

74 

1460. 

0 . 337660E-1 5 

3.16435 

-15.47152 

75 

1480. 

0.287992E-15 

3.17026 

-15.54062 

76 

1500. 

0. 246604E-1 5 

3.17609 

-15.60800 

77 

1520. 

0. 212491 E-15 

3.18184 

-15.67266 

78 

1540. 

0 . 1 84663E -1 5 

3.18752 

-15.73362 

79 

1560* 

0.1 62176E-15 

3.19312 

-15.79001 

80 

1580. 

0.144157E-15 

3.19866 

-15.84116 

81 

1600. 

0.129825E-15 

3.20412 

-15.88664 

82 

1620. 

0.1 1B488E-15 

3.20951 

-15.92632 

83 

1640* 

0.109553E-1S 

3.21464 

-15.96038 

84 

1660. 

0. 10251 7E-15 

3.22011 

-15.98920 

85 

1680. 

0.969663E-16 

3.22531 

-16.01337 

86 

17O0o 

0.925619E-16 

3.23045, 

-16.03357 

87 

1720. 

0.890328E-16 

3.23553 

-16.05045 

88 

1740. 

0.861656E-16 

3.24055 

-16.06465 

89 

1760. 

0.837940E-16 

3.24551 

-16.07678 

90 

1780. 

0.817916E-16 

3.25042 

-16.08728 

91 

1800. 

0 . 800625E-16 

3.25527 

-16.09656 

92 

1820. 

0.785354E-16 

3.26007 

-16.10492 

93 

1840. 

0.771579E-16 

3.26482 

-16.11261 

94 

I860. 

0.758919E-16 

3.26951 

-16.11980 

95 

1880. 

0.747102E-16 

3.27416 

-16.12660 

96 

1900. 

0.735933E-16 

3.27875 

-16.13316 

97 

1920. 

0.725276E-16 

3.28330 

-16.13950 

98 

1940. 

0.715033E-16 

3.28780 

-16.14566 

99 

1960. 

0.705136E-16 

3.29226 

-16.15172 

100 

1980. 

0.695540E-16 

3.29667 

-16.157 67 

101 

2000. 

0.686211E-16 

3.30103 

-16.16353 

102 

2020. 

0.677123E-16 

3*30535 

-16.16933 

403 

2040. 

0.668260E-16 

3.30963 

-16.17505 

104 

2060. 

0.659607E-16 

3.31387 

-16.18071 

105 

2080. 

0. 651 155E-16 

3*31806 

-16.18631 

106 

2100. 

0.642894E-16 

3.32222 

-16.19185 

107 

2120. 

0.634817E-16 

3.32634 

-16.19734 

108 

2140. 

0.626917E-16 

3.33041 

-16.20279 

109 

2160. 

0.619188E-16 

3*33445 

-16.20818 

110 

2180. 

0.61 1625E-16 

3.33846 

-16.21350 

111 

2200. 

0.604223E-16 

3.34242 

-16.21880 

112 

2220. 

0.596975E-16 

3.34635 

-16.22403 

113 

2240. 

0.589679E-16 

3.35025 

-16,22923 

114 

2260. 

0.582929E-16 

3.35411 

-16.23438 

115 

2280. 

0.576121E-16 

3.35793 

-16.23947 

116 

2300. 

0.569451E-16 

3.36173 

-16.24454 

117 

2320. 

0.562915E-16 

3.36549 

-16.24954 

118 

2340. 

0.556509E-16 

3.36922 

-16.25452 

119 

2360. 

0.550230E-16 

3.37291 

-16.25945 

120 

2380. 

0.544074E-16 

3.37658 

-16,26433 

121 

2400, 

0.538037E-16 

3.38021 

-16.26918 

122 

2420. 

0.5321 16E-16 

3.38381 

-16.27399 

123 

2440. 

0.526309E-16 

3.38739 

-16.27875 

124 

2460. 

0.52061 IE-16 

3.39093 

-16.28348 

125 

2480. 

0.515021E-16 

3.39445 

-16.28816 

126 

2500. 

0.509534E-16 

3.39794 

-16.29282 

127 

2520. 

0.504150E-16 

3.40140 

-16.29744 

128 

2540. 

0.498864E-16 

3.40483 

-16.30200 

129 

2560. 

0.493674E-16 

3.40824 

-16.30655 

130 

2580. 

0.488579E-16 

3.41162 

-16.31107 

131 

2600. 

0.483574E-16 

3.41497 

-16.31552 

132 

2620. 

0.478659E-16 

3.41830 

-16.31996 

133 

2640. 

0.473830E-16 

3.42160 

-16.32437 

134 

2660  • 

0.4690B5E-16 

3.42488 

-16.32874 

135 

2680* 

0.464423E-16 

3.42813 

-16.33308 

136 

2700* 

0.459842E-16 

3*43136 

-16.33739 

137 

2720. 

0.455340E-16 

3.43457 

-16.34166 

138 

2740. 

0. 45091 4E- 16 

3.43775 

-16.34590 

139 

2760* 

0.446562E-16 

3.44091 

-16.35011 

140 

2780* 

0..442284E-16 

3.44405 

-16.35429 

141 

2800. 

0.438076E-16 

3.44716 

-16.35844 

142 

2820. 

0.433939E-16 

3.45025 

-16.36256 

143 

2840. 

0.429869E-16 

3.45332 

-16,36665 

144 

2860. 

0.425866E-16 

3.45637 

-16.37073 

145 

2880. 

-0.421927E“16 

3.45939 

-16.37476 

146 

2900. 

0.418052E-16 

3.46240 

-16.37075 

147 

2920. 

0. 414238E-16 

3.46538 

-16,38274 

148 

2940. 

0.410485E-16 

3.46835 

-16.38669 

149 

2960. 

0.40679 IE- 16 

3.47129 

-16.39063 

150 

2980. 

0.403155E-16 

3.47422 

-V 6. 39452 

151 

3000. 

0.399576E-16 

3.47712 

-2.6,39339 

Appendix  B.  Fortran  IV  program  for  C«  model  for  stable 
stratification  with  example  for  2113  MST, 

3  Sept.  1984  for  White  Sands,  NH. 


LIST  CNSQ. STABLE 


1 

£XP0=-4./3. 

2 

WRITE (&» 101 i 

3 

101 

FORMAT (' SEnter  wind  speed  at  10  metres  (in  Knots) 

4 

READC5» 1 02)  WS 

5 

102 

FORMAT <F4 ,1 ) 

6 

WS=US*. 51479 

7 

CNSGC*! » 28£-13*EXP<  <-(WS-4» )4#2)/<2.*2.3»*2) ) 

■8 

WRITE ( 7» 100 )  CNSOC 

9 

100 

FORMAT (E14  *5) 

10 

HGT=10» 

11 

DO  10  1=1 »31 

12 

CN2C=CNSaC*(HGT/4, >**EXPO 

13 

HGTL0G=AL0G10(HGT) 

14 

CN2LOG=ALOG10(CN2C ) 

15 

URI TE  < B » 200 )  HOT »CN2C » HGTLOG » CN2L0G 

14 

200 

FORMAT (F5oOtE15«6»2Xt2F12*5) 

17 

HGT=I*20. 

IS 

10 

CONTINUE 

19 

STOP 

20 

END 

EXAMPLE 


♦RUN  *FTN  SCARDS=CNSQ. STABLE 
♦Execution  begins  15 J 55 l 13 

in  MAIN 

terminated  15.55.14 
7=SURFACE  8=F'R0FILE 
begins  15:55525 
speed  at  10  metres  (in  knots) 


No  errors 
♦Execution 
♦RUN  -LOAD 
♦Execution 
Enter  wind 


J-0,066  $0*04 


♦Execution 

♦ 


terminated  15:55130  T=0,022 


(F4 ♦ 1 ) 

$0.02 


9.9 


2113  MST,  9/3/84 


♦LIST  SURFACE 

1  0.1 1425E-12 


PR0FILE  2 (m) 

cym) 

log1Q(z) 

1o9l0  C* 

1 

10. 

0.336726E-13 

1.00000 

-13.47272 

o 

20. 

0. 133630E-13 

1.30103 

-13.87409 

3 

40. 

0.530312E-14 

1.60206 

-14.27547 

4 

60. 

0.308847E-14 

1.77815 

-14.51026 

5 

80. 

0 . 210454E-14 

1.90309 

-14.67684 

6 

100. 

0 . 156295E-1 4 

2.00000 

-14.80605 

7 

120. 

0. J22566E-14 

2.07918 

-14.91163 

a 

140. 

0 . 997945E- 15 

2.14613 

-15.00089 

9 

160. 

O.835109E-15 

2.20412 

-15.07822 

10 

180. 

0.71  3808E'i5 

2.25527 

-15.14642 

li 

200. 

0.620256E-15 

2.30103 

-15.20743 

12 

220. 

0.546237E-15 

2.34242 

-15,26262 

13 

240. 

0.486403E-15 

2.38021 

-15.31300 

14 

260. 

0 • 437166E-15 

2.41497 

-15.35935 

15 

200. 

0 . 396035E-15 

2.44716 

-15.40227 

1A 

300. 

0.3&1229E-15 

2.47712 

-15.44221 

17 

320. 

0.331445E-15 

2.50515 

-15,47959 

18 

340. 

0 . 305708E-15 

2.53148 

-15.51469 

19 

360. 

0.283275E-15 

2 . 55630 

-15.54779 

20 

380. 

0 . 263572E-15 

2.57978 

-15.57910 

21 

400. 

0.246149E-15 

2.60206 

-15.60880 

w 

420. 

0 . 230646E-15 

2.62325 

-15.63705 

23 

440. 

0.216774E-15 

2.64345 

-15.66399 

24 

460. 

0 . 204300E-15 

2.66276 

-15.68973 

25 

480. 

0. 1°3029E-15 

2,68124 

-15.71438 

26 

500. 

0.182903E-15 

2.69897 

-15.73801 

27 

520. 

0. 173490E-15 

2.71600 

-15.76072 

28 

540. 

0 . 164976E-15 

2.73239 

-15.78258 

29 

560. 

0. 157167E-15 

2.74819 

-15.00364 

30 

580. 

0.1 49982E-15 

2.76343 

-15.82396 

31 

600. 

0. 14335 4E-15 

2 . 778i5 

-15.84359 

